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PREFACE 


This  document  has  been  submitted  under  Contract  Number  DACA72-03-C-0001.  The 
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Summary 


This  report  provides  an  overview  of  the  results  for  the  base  program  under  SERDP  Project 
Number  UX-1326.  The  base  program  for  this  project  has  been  successfully  completed  and  all 
technical  milestones  have  been  met.  The  focus  of  efforts  during  the  base  program  has  been  on 
(1)  adaptation  of  the  instrumentation  and  probe  electronics,  (2)  testing  of  an  enhanced  drive 
winding  method  developed  under  JENTEK  IR&D,  (3)  fabrication  and  testing  of  the  elements  of 
a  sixteen  channel  array  with  higher  resolution  than  the  original  prototype  MWM®- Array,  (4) 
refinement  and  validation  of  JENTEK’ s  model  based  inversion  method  for  estimation  of  depth 
and  size  of  spherical  objects,  (5)  preliminary  field  testing,  (6)  development  of  UXO  depth  and 
size  measurement  method,  and  (7)  preliminary  development  of  model-based  clutter  suppression 
methods. 

Significant  results  achieved  during  the  base  program  have  included: 

•  Demonstration  of  a  substantial  improvement  in  signal  to  noise  and  associated  depth  of 
sensitivity  for  buried  bomblets  and  fabricated  spherical  objects.  This  was  accomplished 
through  an  enhanced  drive  winding  construct  using  adapted  impedance  instrumentation  and 
probe  electronics  developed  under  JENTEK  IR&D. 

•  Demonstration  that  doubling  the  resolution  of  the  sensing  element  array  can  be 
accomplished  without  reducing  signal-to-noise.  Thus,  higher  resolution  images  can  be 
provided  without  loss  of  depth  of  sensitivity.  This  higher  resolution  image  will  be  used  to 
suppress  clutter  and  improve  object  feature  recognition. 

•  Demonstration  of  rapid  scanning  capability.  This  permits  real-time  data  acquisition  with 
images  created  at  scan  rates  of  up  to  4  ft/sec  with  single  frequency  operation.  Higher  data 
acquisition  rates  and  higher  resolution  images  can  be  achieved  with  higher  frequency 
operation  with  instrumentation  enhancements.  However,  for  object  identification,  lower  or 
multiple  frequency  operation,  that  can  lead  to  lower  scan  speeds,  may  be  necessary.  This 
potentially  permits  real-time  analysis  for  classification  based  on  estimated  object  properties. 

•  Demonstration  of  preliminary  object  shape  discrimination  capability.  Preliminary  scan 
images  over  simple  objects,  such  as  a  spherical  bomblet  and  a  cylinder,  indicate  object  shape 
and  orientation  (in  a  plane  parallel  to  the  sensing  array). 

•  Demonstration  of  preliminary  model  based  depth  and  size  estimation  potential  using  a 
model  based  grid  method.  TTiese  methods  provide  estimates  of  object  properties  to  assist 
in  object  identification  and  clutter  suppression. 

•  Demonstration  of  object  imaging  in  the  field.  Preliminary  measurements  were  performed 
with  the  prototype  drive  and  adapted  instrumentation.  Object  responses  obtained  in  the  field 
were  similar  to  laboratory  measurements. 

•  Generated  criteria  for  enhancements  and  recommendations.  Observations  from  the  field 
and  laboratory  measurements  have  been  incorporated  into  enhancement  criteria  and 
recommendations  for  the  follow-on  program. 
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Background 


The  efficiency  of  unexploded  ordnance  (UXO)  remediation  is  currently  limited  by  the  inadequate 
discrimination  capability  of  existing  UXO  detection  technologies.  Inductive  sensors  with  a 
single  sensing  coil  are  commonly  used  for  detection.  These  single  coil  methods  often  detect  all 
relevant  metal  objects  but  generally  cannot  discriminate  between  UXO  and  harmless  clutter.  The 
goal  of  this  program  is  to  develop  inductive  arrays  for  high  resolution  UXO  imaging  and 
discrimination  by  leveraging  existing  MWM- Array  sensor  and  instrumentation  technologies  used 
in  nondestructive  testing. 

Numerous  researchers  have  demonstrated  the  value  of  increased  resolution  through  simulation 
and  field  testing  with  single  coils  scanned  carefully  to  produce  registered  images  of  object  shape. 
JENTEK  has  been  working  on  such  solutions  for  many  years  and  has  fielded  systems  in 
traditional  nondestructive  testing  applications  that  use  a  combination  of  high  resolution  inductive 
sensors  and  model  based  inversion  methods  to  estimate  multiple  unknowns  of  subsurface  cracks, 
and  corrosion  damage. 

JENTEK’ s  MWM- Array  technology  is  based  on  unique  designs  for  EMI  sensor  arrays,  which 
use  a  single  drive  winding  with  multiple  sense  elements,  which  permit  both  high  resolution 
imaging  and  accurate  modeling  of  the  sensor  response  to  buried  objects.  The  drive  creates  a 
shaped  magnetic  field  pattern,  which  provides  a  continuous  variation  in  the  orientation  of  the 
magnetic  field  relative  to  the  buried  ordnance.  Arrays  of  small  inductive  coils  placed  throughout 
the  shaped  field  sense  the  field  variations  from  conducting  or  magnetic  UXO  and  clutter.  Images 
obtained  from  scans  over  buried  objects  provide  a  basis  for  spatial  filtering  and  signal 
processing.  An  eight-element  prototype  array  developed  five  years  ago  demonstrated  the 
capability  to  produce  images  of  buried  ordnance.  Although  this  was  at  a  much  higher  resolution 
than  single  coil  methods,  JENTEK’ s  new  commercial  parallel  architecture  instrumentation  can 
acquire  data  from  up  to  39  channels  in  less  than  0.01  seconds. 

JENTEK  has  also  developed  methods  that  use  pre-computed  databases  of  sensor  responses  to 
analyze  the  response  of  the  sensors  in  real-time.  This  is  critical  since  real  time  imaging  in  the 
field  is  necessary  if  confirmatory  methods  are  to  be  employed  at  the  UXO  site.  If  successful,  this 
program  will  provide  an  enhanced  capability  for  the  reclamation  of  lands  contaminated  with 
unexploded  ordnance.  The  ability  to  discriminate  between  ordnance  and  clutter  and  provide 
information  on  ordinance  type,  depth  and  orientation  would  greatly  reduce  the  false  alarm  rate 
and  improve  the  efficacy  of  reclamation  efforts,  substantially  reducing  the  cost  and  risks 
associated  with  UXO  clearance  and  transitioning  of  lands  back  to  public  use. 
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Results 


This  section  provides  a  summary  of  the  results  from  this  program. 

For  reference,  the  base  program  objectives  were 

(1)  adapt  JENTEK's  commercial  39-channel  impedance  instrument  for  use  with  large  scale 
MWM- Arrays  for  UXO  detection 

(2)  establish  MWM- Array  object  size  and  depth  measurement  capabilities  for  simple  objects 

(3)  evaluate  current  UXO  discrimination  capabilities,  and 

(4)  complete  an  initial  field  test  of  the  prototype  MWM-Array  with  adapted  instrumentation. 

The  tasks  for  the  base  program  plan  were: 

Task  1 :  Adapt  current  instrumentation 

Task  1.1:  Adapt  prototype  array  and  perform  initial  laboratory  tests 
Task  1.2:  Adapt  object  size,  shape,  and  depth  algorithms 
Task  1.3:  Fabricate  12  or  more  sense  element  array 
Task  1.4:  Test  new  array  and  evaluate  performance  in  laboratory  tests 
Task  1.5:  Evaluate  trade-offs  based  on  sense  element  size 
Task  2:  Demonstrate  prototype  array 

Task  2.1:  Acquire  representative  inert  munitions 

Task  2.2:  Develop  initial  detection/suppression  algorithm  for  field  testing 
Task  2.3:  Demonstrate  size  and  depth  measurements  for  known  objects 
Task  2.4:  Demonstrate  preliminary  discrimination  and  clutter  suppression  approach 
Task  2.5:  Provide  plan  for  initial  field  test 
Task  2.6:  Perform  initial  field  tests 
Task  3  :  Determine  enhancement  criteria  and  recommendations 
Task  4:  Submit  Base  Program  final  report 

The  following  sections  include:  (1)  a  description  of  the  system  adaptations;  (2)  laboratory  testing 
results;  (3)  field  testing  results;  and  (4)  a  discussion  of  recommendations  and  criteria  for 
enhancements. 
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System  Adaptations 

As  part  of  this  SERDP  effort,  JENTEK  adapted  the  prototype  instrumentation  to  support  a 
modified  winding  design  that  provides  higher  sensitivity  to  buried  objects.  This  has  included 
adaptation  of  the  MWM- Array  probe  electronics  for  the  prototype  MWM- Array  so  that 
JENTEK’ s  parallel  architecture,  39  channel  impedance  instrumentation  can  be  used  to  support 
rapid  imaging.  Additional  field  hardening  of  the  probe  electronics  prototype  and  MWM- Array 
will  be  required  to  support  extensive  field  testing,  but  the  current  configuration  was  successfully 
used  to  support  laboratory  testing  and  limited  field  testing  during  the  base  program. 

To  support  the  development  of  an  enhanced  imaging  capability,  16  sensing  elements  were 
fabricated.  This  essentially  doubled  the  imaging  resolution  capability  of  the  MWM- Array  over 
the  previous  prototype.  The  array  is  designed  to  fit  within  the  0.61  m  (24  in.)  width  of  the 
existing  drive  winding  structure  with  a  reasonable  spacing  between  the  ends  of  the  array  and  the 
retum/side  portions  of  the  drive  winding.  The  distance  to  the  return  portions  of  the  drive 
winding  needs  to  be  large  to  minimize  extra  coupling  of  the  magnetic  field  between  the  drive 
winding  and  the  sense  coils  since  this  extra  coupling  is  a  source  of  element-to-element  variability 
for  the  array  response.  Figure  1  shows  the  primary  winding  structure  with  the  lower  resolution  1 
in.  square  sensing  elements  configured  in  an  eight-element  array.  Figure  2  shows  one  of  the  new 
sense  elements. 


Figure  1.  Photograph  of  the  “old”  sensing  array  consisting  of  2.54  cm  (1  in.)  square  sense  element  coils  located 
within  the  primary  winding  structure.  The  beam  suspended  above  the  sensor  was  used  for  initial  manual  scans  of 
the  blue  bomblet  (right)  and  a  7.6  cm  (3  in.)  diameter  aluminum  sphere  (left). 
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Figure  2.  Photograph  of  a  “new”  12.7  cm  (0.5  in.)  diameter  sense  element  coil  placed  in  the  drive  winding  structure 
for  initial  testing. 


A  second  linear  array  of  16  sense  coils  is  being  considered  for  fabrication  in  the  follow-on 
program.  The  arrays  would  then  be  placed  at  different  locations  in  the  drive  winding  footprint 
and  could  provide  complementary  information  about  the  object  properties.  The  second  array 
could  also  be  placed  next  to  first  array,  but  offset  a  half-sense  element  in  the  array  direction,  to 
provide  an  even  higher  resolution  image.  Other  sense  coil  orientations  (such  as  vertical  versus 
horizontal)  can  also  be  used  to  provide  complementary  information  about  the  object  properties. 
The  tradeoffs  associated  with  these  various  approaches,  such  as  the  information  gained  from  the 
additional  elements  versus  the  increased  complexity  of  the  electronics,  will  be  explored  further  in 
the  follow-on  program.  Also,  in  the  follow-on  effort  we  expect  to  use  a  second  drive  winding 
oriented  perpendicular  to  the  original  drive  winding  as  conceived  by  Goldfine,  et.al.,  in  the  mid 
1990’s  for  UXO  and  landmine  detection. 


Initial  testing  was  performed  on  a  single  element  of  the  array.  Initial  tests  with  single  sensing 
elements  have  demonstrated  improved  sensitivity.  Three  single  element  configurations  were 
tested:  the  “old”  2.54  cm  (1  in.)  square  coil,  the  “new”  12.7  cm  (0.5  in.)  diameter  round  coil  and 
the  “new”  coil  with  an  enhanced  drive  winding  construct.  A  model  for  estimating  the  signal-to- 
noise  ratio  (SNR)  was  adapted  to  the  MWM- Array  and  indicated  that  modest  modifications  to 
the  drive  winding  can  help  enhance  sensitivity  to  buried  objects.  The  two  objects  used  in  the 
initial  tests,  as  shown  in  Figure  1,  were  a  7.6  cm  (3  in.)  diameter  aluminum  sphere  and  a  6.4  cm 
(2.5  in.)  diameter  blue  bomblet  from  a  cluster  bomb.  These  objects  were  manually  scanned 
along  the  beam  (see  Figure  12)  at  a  rate  of  approximately  0.46  m/s  (1.5  ft/s)  over  the  various 
single  sensing  element  configurations.  Figure  2  shows  the  placement  of  a  “new”  sensing 
element  coil  placed  at  the  center  of  the  drive  winding. 


Measurements  were  made  as  the  bomblet  and  aluminum  sphere  were  scanned  over  the  sensing 
element  at  various  heights.  These  results  are  expressed  in  terms  of  the  signal-to-noise  ratio 
determined  from 


SNR  = 


(  7  _  7 
^ib 

{  AZr  J 

l  AZ,  J 

(1) 


Page  8  of  66 


where  Z  is  the  measured  impedance,  the  subscript  r  denotes  the  real  part,  the  subscript  i  denotes 
the  imaginary  part,  the  subscript  b  denotes  the  baseline  response,  and  A  denotes  the  noise  level. 
The  peak  SNR  value  versus  depth  is  plotted  in  Figure  3,  with  individual  scans  plotted  in  Figure  4 
for  the  bomblet  and  Figure  5  for  the  aluminum  sphere.  In  each  case  the  peak  in  the  SNR  level 
occurs  when  the  object  is  directly  over  the  sense  element.  The  variation  in  the  position  of  the 
peak  is  a  result  of  the  manual  scanning  of  the  object. 

These  plots  and  the  measurement  scan  show  that  the  “old”  2.54  cm  (1  in.)  square  coil  and  the 
“new”  12.7  cm  (0.5  in.)  diameter  round  coil  have  similar  sensitivity  to  the  objects.  This  is 
expected  since  the  spatial  variation  of  the  secondary  fields  from  objects  of  these  sizes  and  at 
these  distances  will  not  be  substantially  different  for  sensing  elements  of  this  size.  This 
indicates  that  there  is  no  signal  degradation  due  to  the  use  of  the  smaller  coils;  the  benefit  of  the 
smaller  coils  will  be  in  the  generation  of  higher  resolution  spatial  images  that  should  improve 
object  discrimination.  The  use  of  the  enhanced  drive  winding  construct  provides  a  substantial 
improvement  in  the  sensitivity  for  both  shallow  and  deep  objects.  The  improvement  in 
sensitivity  is  shown  distinctly  in  several  selected  scans  of  Figure  4a  and  Figure  5a,  where  scans 
with  the  enhanced  drive  winding  construct  gives  a  larger  signal  at  a  deeper  depth  than  the 
previous  drive  winding. 
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Figure  4.  Signal-to-noise  ratio  plots  for  the  blue  bomblet.  The  scan  position  is  the  position  of  the  bomblet  as  it  is 
scanned  over  the  sensing  element  The  dashed  vertical  line  indicates  the  nominal  position  of  the  sense  element 
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Figure  5.  Signal-to-noise  ratio  plots  for  the  3  in.  diameter  aluminum  sphere.  The  scan  position  is  the  position  of  the 
sphere  as  it  is  scanned  over  the  sensing  element  The  dashed  vertical  line  indicates  the  nominal  position  of  the  sense 
element. 
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The  incorporation  of  the  sensor  array  into  a  cart  for  scanning  over  objects  is  shown  in  Figure  6. 
For  rapid  imaging  and  measurements,  individual  connections  are  made  to  each  individual  sense 
coil.  JENTEK’s  commercial  parallel  architecture  impedance  instrumentation  can  accommodate 
the  simultaneous  measurement  of  up  to  39  measurement  channels.  The  coils  are  connected  to  the 
instrumentation  through  an  interface  or  probe  electronics.  These  electronics  provide 
amplification  of  the  measurement  signals  and  some  signal  filtering.  Since  the  measurements 
with  these  sense  coils  and  drive  winding  require  gain  settings  different  from  our  standard  probe 
electronics,  different  internal  circuit  boards  were  used  for  the  probe  electronics.  The  data 
acquisition  is  controlled  via  a  laptop  computer.  In  this  version,  the  laptop  and  support 
instrumentation  were  not  placed  onto  the  cart  structure  itself.  Future  designs  are  expected  to 
have  the  support  electronics  and  laptop  in  a  more  portable  configuration,  either  incorporated  into 
the  cart  itself  or  placed  into  backpacks. 


Figure  6.  Photograph  of  the  sensor  array  mounted  onto  a  cart  for  scanning  over  objects  and  also  the  support 
instrumentation. 
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Models  for  system  response 

Modeling  efforts  have  focused  on  improved  representation  of  the  drive  windings,  the  interactions 
of  multiple  objects,  and  the  effects  of  the  surrounding  media  (soil)  on  the  measurement  response. 
New  models  account  for  both  the  drive  winding  enhancement  and  the  effects  of  the  return 
windings  not  addressed  in  our  previous  model1’2 3.  Accurate  models  of  the  magnetic  fields  created 
by  the  drive  winding  structure  have  been  used  to  calculate  the  secondary  dipole  fields  excited  for 
objects  which  can  be  represented  by  a  polarization  tensor;  spherical  objects  have  been 
specifically  used  in  initial  modeling  and  experiments.  The  sensing  element  response  is  then 
calculated  from  the  secondary  fields  for  spheres  of  various  diameters,  positions  (depth  and 
translation)  and  electrical  properties  (conductivity  and  permeability).  Responses  for  spheres  of 
fixed  electrical  properties  located  directly  over  the  sensing  element  at  varying  depth  and  with 
varying  diameter  have  been  tabulated  into  depth-diameter  measurement  grids,  allowing  quick 
inversion  from  measured  sensing  element  response  (transimpedance)  to  depth-diameter 
parameters. 

These  models  extend  previous  work  by  permitting  interactions  between  multiple  objects 
(spheres)  and  by  permitting  more  realistic  field  calculations  for  the  drive  winding.  The  modeling 
approach  here  follows  the  approach  of  Braunisch  et  at  for  distributions  of  conducting  and 
permeable  bodies  of  generic  geometry  in  the  presence  of  an  imposed  magnetoquasistatic  field. 
The  model  accounts  for  the  interactions  of  secondary  fields  generated  by  the  objects.  Solving 
this  multi-region,  arbitrary  geometry  problem  with  an  exact  analytic  solution  is  currently 
impractical.  Therefore  the  problem  is  simplified  in  several  ways  such  that  an  approximate 
solution  can  be  obtained. 

The  system  is  assumed  to  be  magnetoquasistatic,  which  amounts  to  ignoring  the  wave  nature  of 
the  fields  by  allowing  the  displacement  current  to  be  ignored  in  Maxwell’s  equations.  The 
conductivity  of  the  material  separating  the  objects  is  assumed  to  be  weakly  conducting  such  that 
only  the  objects  induce  secondary  magnetic  fields.  These  two  requirements  can  be  expressed 
mathematically  by  requiring  the  magnitude  of  the  complex  wavelength  to  be  much  larger  than 
the  length  scales  of  the  system.  The  length  scale  of  the  system  includes  sizes  of  objects,  distance 
to  the  sources,  distance  to  observation  points,  and  inter-object  distance.  The  weakly  conducting 
criteria  also  require  the  conductivity  of  the  objects  to  be  significantly  larger  than  the  background 
medium  when  the  object  has  low  permeability.  Alternatively  the  permeability  of  the  object 
should  be  significant  for  weakly  conducting  objects.  Though  difficult  to  quantify,  this  criteria 
helps  assure  that  the  object  effect  will  be  dominant  over  the  diffusion  effects. 

The  fields  are  also  assumed  to  be  spatially  uniform  over  the  region  of  each  object.  This  requires 
the  variation  of  the  fields  due  to  both  the  primary  sources  and  due  to  the  induced  fields  from 
neighboring  objects,  to  be  reasonably  constant  over  the  region  of  each  object.  The  dominant 


1  N.  Goldfine,  A.  Washabaugh,  D.  Schlicker,  “High-Resolution  Inductive  Array  Imaging  of  Buried  Objects,”  SPIE 
AeroSense  Conference,  Orlando,  FL;  April  2000. 

2  N.  Goldfine,  D.  Schlicker,  A.  Washabaugh,  D.  Clark,  M.  Zahn,  “New  Quasistatic  Magnetic  and  Electric  Field 
Imaging  Arrays  and  Algorithms  for  Object  Detection,  Identification,  and  Discrimination,”  SPIE  AeroSense 
Conference,  Orlando,  FL;  April  1999.  Proceedings  vol.  3710;  n.1, 1999,  pp.  89-100. 

3  Henning  Braunisch,  Chi  O.  Ao,  Kevin  O’Neill,  and  Jin  A.  Kong,  “Magnetoquasistatic  Response  of  a  Distibution  of 
Small  Conducting  and  Permeable  Objects,”  Geoscience  and  Remote  Sensing  Symposium,  2000.  Proceedings,  vol.  4, 
2000,  pp.  1424-1426. 
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field  from  each  object  is  assumed  to  be  that  of  the  far-field  of  a  dipole.  Most  arbitrary  objects 
would  require  higher  order  multipole  responses  for  an  exact  representation,  even  with  uniform 
excitation.  Fortunately,  these  higher  order  poles  have  increasingly  faster  spatial  decays  and 
therefore  can  be  neglected  at  sufficient  distance.  The  dipole  response  has  a  derivative,  which 
decays  with  distance  and  therefore  becomes  more  spatially  uniform.  This  allows  the  uniformity 
assumption  to  be  achieved  by  forcing  object  spacing  to  be  sufficiently  large  such  that  all  higher 
pole  responses  are  negligible  and  the  dipole  spatial  dependence  has  slowed.  Neglecting  the 
higher  order  poles  also  requires  that  the  fields  only  be  observed  at  sufficient  distance  from  the 
object,  for  objects  which  do  not  have  a  pure  dipole  response. 


The  sources  induced  in  the  objects  are  in  general  made  up  of  both  magnetization  and  diffusion 
currents.  These  currents  are  not  necessarily  in  phase  with  the  excitation  field  and  therefore  the 
dipole  moment  approximating  these  sources  will  be  complex.  The  far  field  response  of  a  z 
directed  dipole  located  at  the  origin  in  spherical  coordinates  is  as  follows: 

H  =  -^L-(2cos6,ir+sin#if?)  (2) 

4  nr 

By  using  the  following  Spherical-Cartesian  coordinate  relations: 
x  =  rsin#cos^ 
y  =  rsin#sin^ 

z-r  cos#  P) 


ir  =  sin  6  cos  +  sin  6  sin  <jAy  +  cos^i* 
L  =  cos#cos0iz  +  cos# sin -sin#iz 


Eq.  (2)  can  be  expressed  as: 


H= 


4  nr 


Jxz.  y  z .  z  z . 

3 - 1  H - i  h - iz 

yr  r  r  r  r  r  ) 


(4) 


The  dipole  moment  of  an  object  may  have  any  orientation,  subject  to  imposed  fields  and 
additional  relations,  which  can  be  broken  down  into  its  mutually  orthogonal  components.  The 
fields  from  the  two  additional  components  can  be  obtained  simply  by  applying  the  following 
coordinate  rotations  to  Eq.  (4): 

z  —>x,it  — » ix  z~*y, 

y  ->  y,  iy  ->  iy  x  ->  x,  ix  ->  ix  (5) 

x ^  i  — ^  — iz  y  — ^  2,  iy  — ^  — iz 

The  location  of  the  dipole  can  be  shifted  to  a  location  rn  =  x„iz  +^niy  +  z„i,  while  the  fields  are 


observed  at  rm  =  xmiz  +^miy  +  zmiz  so  that: 

x=xm-xn,y  =  ym-yH>z  =  zm-zn  (6) 

After  applying  the  coordinate  rotations,  shifts,  and  summing  the  fields  from  each  component  of 
the  dipole  moment,  the  total  H-field  can  be  decomposed  into  components  as  follows: 


4  nr 


v  r  r  / 


-+3m. 

r  r  r  r 


(7) 
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where 


H, 


1 


y  W 
1 


H  = 


!  W 


3m  (y"  y")  (x-  i  m  y»)(-y”  ^)_il  +  3m  y")^Z"  *») 

1  r  r  y\  r  r  J  r  r 

3m  (z”  ~Z”)  (*"  ~X")  |  3m  (z"  ~Z")  0^ ~^»)  j  m  3(^2 _ Z")  (Z" — 5il_] 

1  r  r  y  r  r  z{  r  r 


(8) 

(9) 


which  is  written  compactly  as: 

3(im„ )T  i„ 


H„(r.)  = 


4 jtri 


where 


(10) 

(11) 


1  m„  ;  ^=km-^|  (12) 

Eq.  (11)  is  utilized  in  determining  the  contribution  to  the  fields  at  object  m  of  object  n  based  on 
object/?  ’s  effective  dipole  moment  mn . 


The  effective  dipole  moment  of  an  object  is  dependent  on  the  magnetization  and  diffusion 
currents  present  in  the  object,  which  are  due  to  the  external  excitation  fields.  Since  these 
excitation  fields  have  been  assumed  uniform  over  the  object’s  volume  and  the  resulting  fields 
from  each  object  is  described  by  a  net  dipole  moment,  a  polarizability  tensor  can  be  used  to 
relate  the  fields  to  the  dipole  moment.  The  generic  polarizability  tensor  has  the  form  below 

where  the  elements  of  the  tensor  M  may  be  complex  due  to  the  phase  of  the  object’s  induced 
current.  Since  diffusion  effects  will  be  frequency  dependent  and  the  objects  permeability  may  be 
dispersive,  the  tensor  may  also  be  a  function  of  frequency. 


Miy 

Hx 

my 

= 

My, 

My. 

Hy 

M„ 

Miy 

M..J 

LH. 

Each  column  of  M  represents  the  dipole  moment  that  would  be  excited  by  uniform  magnetic 
field  in  the  corresponding  coordinate  vector  direction.  Performing  the  matrix  multiplication  on 
the  right  essential  superimposes  the  individual  moments  excited  by  each  vector  component  of  the 


field. 


The  polarizability  tensor  for  an  object  can  be  determined  in  several  ways.  Some  geometries  such 
as  spheres  and  spheroids  allow  for  analytic  solution.  Other  more  complex  geometries  may 
utilize  numerical  techniques  such  as  MAS  (Method  of  Auxiliary  Sources),  TSA  (Thin  Skin 
Approximation),  BEM  (boundary  integral  method),  or  FEM.  Instruments  can  also  be  setup  to 
determine  the  polarizability  tensor  empirically,  allowing  characteristics  of  objects  which  may  be 
difficult  to  simulate,  to  be  captured.  All  these  techniques  allow  a  library  of  tensors  to  be 
compiled,  which  can  be  quickly  accessed  even  if  the  initial  computation  or  measurement 
required  significant  time. 

The  relations  for  the  fields  induced  by  the  objects’  dipole  moments  and  the  relation  for  the 
objects’  dipole  moments  induced  by  the  fields  can  now  be  applied  for  each  object  within  the 
system.  The  primary  magnetic  fields,  from  the  sources  external  to  the  objects,  needs  to  be 
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specified  at  each  object  location.  The  only  requirement  is  that  again  these  fields  are 
approximately  uniform  over  the  region  of  each  object.  It  is  also  assumed  that  the  secondary 
object  fields  will  not  affect  the  sources  of  the  primary  field. 


At  each  object  location  the  polarization  tensor  relates  the  total  magnetic  field  at  the  objects 
location  to  the  object’s  induced  dipole  moment  by  Eq.  (13).  The  total  field  at  an  object’s 
location  is  the  sum  of  the  contributions  from  the  primary  field  and  all  other  objects.  In  block 


matrix  notation  the  equation  is: 


mi 

mN 


M,  0 

o 

0  0 


0 

0 

mn 


o  h,(r,)  •••  hN(r,) 

h,(A 

Si(r2)  0  *\  : 

+ 

•  ’•  0  ^N^AT-l)) 

mN 

Si  (rN)  -  Va)ft)  0 

(14) 


where,  Hj  (r„)is  the  incident  (primary)  field  at  location  rn  of  the  nth  object,  mnis  the  dipole 
moment  of  the  nth  object,  Mn  is  the  polarization  tensor  of  the  nth  object,  and 

(15) 

is  a  tensor  relating  the  magnetic  field  at  location  rm  to  a  dipole  at  rn .  The  zeros  in  the  matrix, 

containing  the  polarization  tensors,  are  present  because  the  dipole  moment  for  each  object  is  not 
directly  dependent  on  the  fields  at  other  object  locations.  The  zeros  in  the  matrix  containing  the 

hn  (rm)  are  present  because  the  induced  dipole  moment  for  each  object  only  contributes  to  the 

field  for  other  objects  and  not  the  object  itself.  The  solution  is  now  obtained  by  rearranging  Eq. 
(14)  into  the  standard  Ax-b  form  of  Eq.  (15)  and  solving  for  the  dipole  moments  of  the  objects. 

1  -Mjh2  (/[)  — 

-mAOs)  i  ^ 

■-  i 

-MNhj  (rN )  -MNh(N1^  (rN )  1 

For  computational  purposes  it  is  useful  to  note  that  hn  (rm)  =  hm  (rn) .  It  is  also  useful  to  examine 

the  affects  of  changing  system  parameters  on  the  system  equation,  for  the  purpose  of  optimizing 
computation  during  parametric  sweeps.  The  following  table  summarizes  the  quantities  effected 
by  changing  the  listed  properties  and  whether  the  A  matrix  and/or  b  vector  of  Eq.  (19),  in  Ax  =  b 
form,  is  affected  and  to  what  extent. 


m. 


m. 


^nH«  (av) 


(16) 


4  nr 


In  all  cases  the  “b  vector”  is  affected,  while  others  only  require  a  change  to  the  “A  matrix”  and  in 
many  cases  only  part  of  the  elements  will  need  to  be  recalculated.  Even  in  the  cases  where  all 
elements  need  recalculation,  many  of  the  components  of  the  products  that  appear  in  the  elements 
of  Eq.  (16),  can  be  stored  to  reduce  repetitive  computation  during  parameter  sweeps. 
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Once  the  dipole  moments  have  been  solved,  it  is  a  simple  matter  to  determine  the  total  field  at 
any  location  r  .  This  is  accomplished  by  applying  Eq.  (1 1)  with  rm  =  r  for  each  object  and 
summing  the  resulting  fields  with  the  primary  field  at  r  . 


Parameter  Changed 

Affected  Quantities 

A  Matrix 
Affected 

b  Vector 
Affected 

nth  Object  Location 

K{rm),  Ht(r„) 

Part 

Part 

nth  Object  Property 

mn 

Part 

Part 

Primary  Source 

Change 

None 

All 

Equal  Translation  of 

All  Objects 

a,  ft) 

None 

All 

Frequency 

All** 

All 

*  H,  (/;)  is  only  changed  if  the  primary  source  has  a  frequency  dependence. 
**  Except  for  the  obviously  constant  diagonal  elements 


Several  simulations  were  performed  to  demonstrate  the  utility  of  the  model  and  the  effects  of 
object  interaction.  In  these  simulations  a  square  coil  was  chosen  as  a  realistic  and  common 
excitation  source.  (A  distributed  drive  winding  has  also  been  simulated  with  this  method  for 
generation  of  the  measurement  grids  below.)  Assuming  that  the  wires  are  infinitely  thin  "current 
sticks"  and  carry  a  constant  current,  then  the  field  can  be  determined  from  the  Biot-Savart  Law 


H  ‘ 


.A 


cxa 

ac 

ab 

1-  — 12 
|cxa| 

|c| 

V  1  1 

b 

7 

(20) 


where,  I  is  the  current  in  the  current  stick,  b  is  a  vector  from  the  observation  point  to  the  start  of 
the  stick,  c  is  a  vector  from  the  observation  to  the  end  of  the  stick  and  a  =  c  -  b .  The  excitation 
field  at  each  object  location  is  obtained  by  summing  the  contribution  from  each  of  the  four  legs 
of  the  square  coil. 


Figure  7  shows  the  resulting  induced  or  secondary  magnetic  field  intensity  parallel  to  the  coil 
axis  and  at  the  center  of  the  coil  as  the  coil  is  scanned  over  one  or  more  spheres.  The  excitation 
frequency  was  15  kHz,  the  coil  was  a  square  with  10  cm  sides,  and  the  sphere  conductivity  was 
107  S/m.  The  configurations  considered  were  1)  a  single  1  cm  diameter  sphere,  2)  a  single  1.48 
cm  diameter  sphere,  and  3)  three  1  cm  diameter  spheres  stacked  vertically  with  a  3  cm  center-to- 
center  distance.  The  effects  of  ignoring  the  interactions  are  especially  noticeable  in  the 
imaginary  part  of  the  response.  The  size  of  the  large  sphere  was  choosen  to  produce  a  real 
response  equal  to  that  of  the  three  interacting  smaller  spheres.  The  imaginary  response  is  very 
different  though  between  these  two  configurations,  indicating  that  there  is  information  that 
can  be  used  to  discriminate  between  a  large  object  and  multiple  smaller  objects  (e.g.,  using 
a  Grid  Method  or  other  model  based  approach). 


4  H.A.  Haus  and  J.R  Melcher,  “Electromagnetic  Fields  and  Energy,”  Prentice  Hall,  Englewood  Cliffs,  NJ,  1989. 
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Figure  7.  Real  and  imaginary  induced  magnetic  field  for  several  sphere  configurations  translated  across  a  square 
coil  operating  at  15  kHz. 


Another  model  has  been  explored  as  a  method  for  determining  the  secondary  magnetic  fields 
produced  by  a  sphere  excited  by  a  locally  uniform  time  harmonic  magnetic  field.  This  model 
follows  the  approach  of  Wait5  and  is  not  limited  to  the  purely  magnetoquasistatic  situation.  The 
restriction  on  the  conductivity  of  the  region  external  to  the  sphere  is  also  less  severe,  which 
allows  some  degree  of  magnetic  diffusion  in  the  external  region  to  be  accounted  for.  The 
restriction  placed  on  the  conductivity,  frequency,  and  permeability  of  the  external  region  is 
completely  based  on  making  the  assumption  of  a  locally  uniform  excitation  field  valid.  This 
means  that  the  spatial  variation  of  the  excitation  due  to  both  the  wave  nature  of  the  fields  and  the 
diffusion  affects  needs  to  be  small  relative  to  the  sphere.  This  is  quantified  in  terms  of  the 
propagation  constant  and  radius  as  \y2R\ «:  1  where  R  is  the  radius  of  the  sphere  and  the 
propagation  constant  is  defined  as 

r  =  {ia/J6)-sfid)  )  ^ 

The  subscript  2  refers  to  the  quantities  in  the  external  region,  while  a  subscript  1  refers  the 
quantities  inside  the  sphere.  There  is  no  restriction  on  the  propagation  constant  on  the  sphere’s 
interior  since  both  the  wave  nature  and  the  diffusion  effects  are  fully  accounted  for  in  the 
solution. 

To  solve  for  the  field  distribution,  the  magnetic  field  is  defined  in  terms  of  a  vector  potential 
Fas: 

H  =  ~(o+  ia>s)F  +  —  V  ( V  •  F)  (22) 

ipuo  v  7 


5  J.R.  Wait,  “A  Conducting  Sphere  in  a  Time  Vaiying  Magnetic  Field,”  Geophysics,  vol.  16,  no.  1,  January  1951, 
pp.  666-672. 
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which  differs  from  the  more  common  vector  potential  representation  of  the  magnetic  flux  density 
B  =  V  x  A.  To  better  understand  Eq.  (22),  consider  Ampere’s  Law  for  time-harmonic  fields  in  a 
source-free  region  of  uniform  anisotropic  conductivity: 

VxH  =J  +  ( icos)E  =  (a  +  icos)E  (23) 

Taking  the  curl  of  Eq.  22  results  in 

VxH  -  -(cr  +  ico£)V x F  (24) 

By  comparison  with  Ampere’s  Law  in  Eq.  (23)  it  can  be  seen  that  the  vector  potential  F  is 
related  to  the  electric  field  by 

E  =  -VxF  (25) 

Any  arbitrary  solenoidal  electric  field  can  be  described  as  the  curl  of  another  vector  field  in  this 
way.  Since  the  regions  of  interest  are  composed  of  uniform  electrical  properties  and  Eq.  (25) 
guarantees  that  the  electric  field  will  be  solenoidal;  the  requirement  that  there  is  neither  volume 
charge  nor  any  divergence  in  the  current  density  will  always  be  satisfied.  From  Eq.  (24),  H  and 
F  can  only  be  equated  to  each  other  within  the  gradient  of  a  scalar  potential  and  a  constant, 
which  has  been  set  to  zero,  as: 

H  =  -(cr + icos)F +VU  (26) 

The  divergence  of  H  must  be  zero  in  the  region  of  constant  magnetic  permeability,  while  the 
divergence  of  F  has  yet  to  be  defined.  With  E  and  H  related  to  F  and  U  by  Eq.  (25)  and  Eq. 
(26),  Ampere’s  Law  is  always  satisfied.  Faraday’s  Law  for  time-harmonic  fields  is  as  follows: 

Vx£  =  -io/iH  (27) 

Substitution  into  Faraday’s  Law  and  the  application  of  a  vector  identity  for  VxVxF  results  in: 

V2F-V(VF)  =  y2F-iconVU  (28) 

Since  the  divergence  of  F  can  be  defined  in  terms  of  the  scalar  U  as: 

U  =  J— V-F  (29) 

IfiO) 

Eq.  (28)  simplifies  to 

V2F  =  y2F  (30) 

which  is  a  single  differential  equation  for  the  vector  F  in  each  region.  However,  before 
proceeding  with  the  solution,  the  problem  can  be  further  simplified  by  analyzing  the  required 
components  of  F  based  on  the  symmetry  of  the  problem  and  boundary  conditions  that  must  be 
met. 

The  imposed  field  is  assumed  as  uniform  and  z  directed  and  therefore  contains  only  f  and  6 
components.  Since  the  fields  are  rotationally  symmetric  about  the  z  axis  all  derivatives  with 

A 

respect  to  4>  will  be  zero.  Therefore  by  using  Ampere’s  Law  in  spherical  coordinates  it  can  be 
shown  that  it  is  necessary  for  the  electric  field  to  contain  a  <j)  component  in  order  to  satisfy 
normal  and  tangential  boundary  conditions.  Faraday’s  Law  can  also  be  applied  to  an  electric 

/A 

field  with  only  a  <t>  component,  which  in  combination  with  the  axial  symmetry,  produces  a 
magnetic  field  with  only  r  and  0  components.  It  is  therefore  reasonable  to  assume  that  a  valid 
solution  can  be  found  with  a  ^  directed  electric  field.  By  analyzing  Eq.  (25),  it  can  be  shown 
that  the  vector  F  must  contain  an  f  or  6  component  in  spherical  coordinates,  or  a  p  or  z 
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component  in  cylindrical  coordinates,  in  order  for  the  associated  electric  field  to  contain  a  (f> 
component. 


For  a  solution,  it  is  convenient  to  choose  an  F  vector  with  a  z  component,  even  though  it  may 
seem  more  preferable  to  choose  an  F  vector  with  spherical  components  rather  than  cylindrical 
components.  However,  the  choice  of  spherical  components  can  be  shown  to  make  the  problem 
more  complex  since  the  vector  Laplacian  of  (30)  in  spherical  coordinates  results  in  multiple 
vector  components.  This  is  true  even  in  the  case  where  F  contains  only  a  single  spherical 
component.  Solving  the  problem  completely  in  cylindrical  coordinates  would  also  prove  to  be 
difficult  since  the  surface  on  which  the  boundary  conditions  are  to  be  imposed  is  described  more 
simply  in  spherical  coordinates  as  is  the  normal  and  tangential  components  of  the  field  on  the 
boundary.  Consequently  both  spherical  and  cylindrical  coordinate  systems  are  used  by  assuming 
an  F  vector  solution  of  the  form 

F  =  F,(r,0)z  (31) 

This  expression  can  now  be  substituted  into  equation  Eq.  (30)  and  the  vector  operators  applied  to 
produce 


1JL 

P  dp 


(  dFA 

l  dp) 


i  ?F, 

dz 


(32) 


Since  the  partial  derivatives  are  with  respect  to  cylindrical  coordinates,  but  Fz  is  a  function  of 
spherical  coordinates,  the  chain  rule  must  be  applied  using  the  following  partial  derivative 
relations  between  cylindrical  and  spherical  coordinate  systems 


dr  dr  d0 

—  =  sin#,  —  =  cos#,  — 
dp  dz  dp 


cos  9 


80 

dz 


sin# 


This  results  in 


2d2Fz  0  dFz  2  2_ 

-2 -+2r  — L  -  r  y  Fz  = 


dO2  dd 


(33) 


(34) 


dr2  dr 

This  can  be  converted  into  two  regular  differential  equations  by  separation  of  variables,  with 
F,(r,0)  =  F,(r)F,(  cosff)  (35) 

The  first  differential  equation,  which  has  solutions  related  to  Bessel  functions,  is 

r2F*  (r)  +  2rF’  (r)  -  (r2y2  +n(n+ l))Fr  (r)  =  0  (36) 

The  second,  which  is  Legendre’s  Equation,  is 

(l  -  x2)F;  (x)  -  2  xFe’  (x) +n(n+l)Fs  (x)  =  0  (37) 

where  x  =  cos# .  The  solutions  to  these  two  differential  equations,  in  the  form  of  product 
solutions  for  F,  are 


F,  M) 


iFrlp.i  cose) 


K„(yr ) 


(38) 


P»(  cos#) 


with  I„  (yr)  and  K„  (yr)  defined  in  terms  of  the  standard  Bessel  Functions  as 


Page  20  of  66 


I„(rr)  =  J*rri21„drr) 

Kn{yr)  =  j2yr/xKn+L(yr) 


In  the  vicinity  of  the  sphere  the  excitation  has  been  assumed  as  a  uniform  z  directed  magnetic 
field  of  magnitude  H0 .  This  can  be  related  to  the  F  vector  using  Eq.  (22)  as 


:  F  Z  —  ' 
1  0 z 


7 - z 

(<t2  +ms2) 


The  subscript  2  is  used  on  the  electrical  properties  to  indicate  the  region  exterior  to  the  sphere. 
Therefore  this  is  the  representation  of  the  locally  uniform  excitation  fields  by  the  F  vector,  valid 
in  the  exterior  region. 

The  total  F  vector  on  the  interior  of  the  sphere  can  now  be  expressed  as  a  superposition  of  the 
solutions  to  Eq.  (32),  which  have  no  singularities  at  the  origin  and  weighted  with  yet 
undetermined  coefficients  bn  as 

Fx  =  F,(r,0)z  =  ^b„^^P„(cos0)iu&z  (42) 

n= 0  r 

where  the  subscript  1  indicates  parameters  on  the  interior  of  the  sphere.  The  total  F  vector  on 
the  exterior  of  the  sphere  can  be  expressed  as  the  superposition  of  the  solutions  to  Eq.  (32), 
which  have  no  singularities  at  infinity  and  weighted  with  yet  undetermined  coefficients  an ,  and 

the  excitation  F  vector  of  Eq.  (41)  as 

F2=F2{r,e)z=  ~~  P„ (cos 0)in2(o  z— — \  H0i  (43) 

r  J  (a2+ia)e2) 

In  order  to  determine  the  coefficients  an  and  bn ,  Eqs.  (42)  and  (43)  must  satisfy  the  boundary 

condition  on  the  surface  of  the  sphere.  The  boundary  conditions  require  continuity  of  normal 
magnetic  flux  density  and  continuity  of  the  tangential  magnetic  field  expressed  respectively  as 

0  (44) 

rx{H,-H2)  1^=0  (45) 

In  order  to  apply  the  boundary  conditions,  the  magnetic  field  must  be  determined  from  Eq.  (22) 
in  terms  of  Fz .  Since  the  spherical  components  of  the  magnetic  fields  are  desired  in  the 

evaluation  of  Eqs.  (44)  and  (45),  F  is  first  represented  in  spherical  components  as 

— — 

F  =  FZ  cos  Or -Fz  sin  6*0  after  which  Eq.  (22)  can  be  evaluated  as 


H  =  -y2  cos 9FZ  +  cos 0 


d2F,  sin#  d2F  sin#  dFr 


drdG 


d6  I  i^co 


2  .  sin#  SF  cos#  #  F  sin#  d  F_  1 

+  y  sin#F, - z-+ - 2 - = - f  - 1 

r  dr  r  drdQ  r2  SO  i/uco 
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On  the  surface  of  the  sphere  at  r  -  R  the  source  fields  have  an  r  component  with  a  cos# 
dependence  and  a  6  component  with  a  sin#  dependence.  It  can  be  shown  by  substituting  either 
Eq.  (42)  or  (43)  into  Eq.  (46),  that  only  the  n  =  0  terms  of  the  series  have  the  correct  # 
dependence  at  r  =  R ,  and  are  the  only  terms  required  for  the  solution.  Since  P0  (cos#)  =  1 ,  Fz 

has  no  #  dependence  and  all  of  the  terms  in  Eq.  (46)  that  contain  partial  derivative  with  respect 
to  #  can  be  removed.  By  utilizing  Eqs.  (44),  (45),  and  (46)  the  boundary  condition  can  be 
expressed  as 


~d2Fx 

1  ns 

2rr 

~d2F2 

[dr2  /,MJ 

_  dr2 

1 

|  1  ^  ylF 

!_[ 

'1  dF2 

F\ 

~  n  M 

r  or  j 

Hi 

r  dr 

-r22F2 


—rlF2 


(47) 


lr=R 


The  coefficients  a0  and  b0  can  now  be  determined  by  simultaneously  solving  the  equations 
produced  by  substituting  Eqs.  (42)  and  (43)  into  Eq.  (47)  and  eliminating  terms  containing 
positive  powers  of  (j^i?) ,  which  are  assumed  to  be  much  less  than  one.  The  resulting 

coefficients  are 

^a2r0(a)-(2^  +  {i2)ai'0(a)+(2H  +//2)/0  (or )-(//, -y2)a2i0(a) 
//,or2/0 (a)  — 2(^/,  — //2)q'/0 (<af)  + 2(//j  — //2)/0(a)— (//j  +2//2)a  /0(«) 

3/^2 


ao  = 


R2Hn 


K0{r2R) 


(48) 


b0=R*H0 


(49) 


fj,a2iZ{a)- 2(m 2(H\  -/r2)/0(a)-(//,  +2//2)«270(a)J 

where  a  =  yxR .  The  functions  /0  and  K0  can  be  represented  in  terms  of  hyperbolic  and 
exponential  functions6  as 

/„(z)  =  sinh(z) 

K„(z)  =  e- 

These  can  be  used  to  further  simplify  the  coefficients  to: 

2 //,  (sinh  a -a  cosh  a )  +  jj2  (sinh  a  -  a  cosh  a  +  a2  sinh  a ) 

2//,  (sinhor-orcosha)-2//2  (sinh  a  -a  cosh  a + or2  sinh  a) 

_ 3^2 _ 

2/*,  (sinh  a -a  cosh  a)-  2//2  (sinh  a  -  a  cosh  a  +  a1  sinh  a ) 

The  magnetic  fields  on  both  the  interior  and  exterior  of  the  sphere  can  be  evaluated  by  utilizing 
relations  (50)  and  (51)  in  Eqs.  (42)  and  (43)  for  Fx  and  F2,  independently  substituting  these 
expressions  into  Eq.  (46),  and  further  expressing  them  in  cylindrical  coordinates  to  produce: 


ao  -  HqR3 


b0  =  H0R 3 


(50) 

(51) 

(52) 

(53) 


6  Abramowitz,  Milton,  and  Irene  A.  Stegan,  Handbook  of  Mathematical  Functions,  Dover  Publications,  1964,  pp. 
443-444. 
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H]  =  b0  -^[sinh  A  (/?,2  +  3)  -  3/?,  cosh  A]/? 


1  3z 


r  K  r 


P  Pl. 


- 1  (sinh  A  -  A  cosh  A)  — —  Pi  sinh  Ai  z 


i  (A+i)- 


5!=a0e  ft  -^(A2+3A+3)p+^ft^-  fe 


z  (55) 


where  A  =  yxr  and  p2  =  y2r  .  The  single  j>  component  of  the  electric  field  can  also  be 
evaluated  from  Fx  and  F2  by  using  Eq.  25  and  expressing  the  results  in  cylindrical  coordinates 
as  follows 

A*  =  bjnyo  4[  A  cosh  A  -  sinh  A  ]  <56) 


E2<p  =  -W  [A  + 1]  (57> 

This  model  can  then  be  used  to  plot  the  magnetic  field  lines  or  induced  current  densities  resulting 
from  the  simulated  response  of  a  sphere.  The  following  simulations  assumed  a  sphere  of  radius 
5  cm,  conductivity  105  S/m,  relative  permeability  1,  and  relative  permittivity  1.  They  also 
assumed  an  external  medium  of  conductivity  0,  relative  permeability  1,  and  relative  permittivity 
1  The  frequency  was  assumed  to  be  15  kHz,  and  the  primary  field  amplitude  was  assumed  to  be 
1  A/m.  This  results  in  a  skin  depth  which  is  approximately  one  quarter  of  the  sphere’s  radius.  It 
should  also  be  noted  that  even  though  the  displacement  currents  were  not  explicitly  set  to  zero 
(i.e.,  sx  =  s2  =  0 ),  their  affect  on  the  fields  and  the  currents  in  the  region  plotted,  is  negligible. 


The  plots  in  Figure  8  through  Figure  1 1  show  the  magnetic  field  lines  for  the  total  field  (primary 
and  secondary)  and  the  internal  current  densities  for  a  cross-section  of  the  sphere.  Each  set  of 
plots  shows  the  field  lines  and  current  distribution  for  a  particular  phase  of  the  periodic  cycle  in 
increments  of  tt/4  radians.  The  plots  only  include  the  phases:  cot  =  0,^/4,  zr/2, 3^/4 ,  however 
the  cycle  is  completed  by  simply  reversing  the  field  and  current  directions  of  the  first  four 
phases. 

The  general  cyclic  behavior  of  the  current  can  be  visualized  as  a  movement  of  a  concentration  of 
(p  directed  current,  originating  at  the  surface,  into  the  sphere.  These  surface  currents  tend  to 
form  a  pocket  of  current  as  they  move  inward  toward  the  sphere’s  center,  while  simultaneously 
becoming  more  diffuse.  The  magnetic  fields  tend  to  circulate  around  these  pockets  of  current 
and  also  move  inward  until  the  centers  of  the  circulations  connect  and  the  cycle  repeats. 

Figure  10  {cot  -  jv/2)  has  a  magnetic  field  distribution  similar  to  that  associated  with  a  dipole 
field,  which  appears  significantly  different  from  the  other  phases  plotted.  This  is  a  result  of  the 
primary  field  intensity  being  zero  for  cot  =  zr/2,3zr/2 ,  while  the  secondary  fields  are  generated 
from  the  circulating  eddy  currents  which  would  appear  as  the  fields  from  a  loop  or  dipole. 
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(a)  (b) 


Figure  8.  Simulated  magnetic  field  lines  (a)  and  current  density  (b)  for  a  nonmagnetic  conducting  sphere  in  a  locally 
uniform  imposed  magnetic  field  at  a  fixed  phase  (cot=0)  of  the  periodic  excitation  cycle. 


Figure  9.  Simulated  magnetic  field  lines  (a)  and  current  density  (b)  for  a  nonmagnetic  conducting  sphere  in  a  locally 
uniform  imposed  magnetic  field  at  a  fixed  phase  (®t=7t/4)  of  the  periodic  excitation  cycle. 
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Figure  10.  Simulated  magnetic  field  lines  (a)  and  current  density  (b)  for  a  nonmagnetic  conducting  sphere  in  a 
locally  uniform  imposed  magnetic  field  at  a  fixed  phase  (cot=7t/2)  of  the  periodic  excitation  cycle. 
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Figure  1 1.  Simulated  magnetic  field  lines  (a)  and  current  density  (b)  for  a  nonmagnetic  conducting  sphere  in  a 
locally  uniform  imposed  magnetic  field  at  a  fixed  phase  (<»t=37i/4)  of  the  periodic  excitation  cycle. 
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Laboratory  Tests 

Numerous  tests  were  performed  in  a  laboratory  setting  to  demonstrate  the  basic  operation  of  the 
sensor  and  grid  measurement  methods.  Initial  validation  of  the  models  for  the  sensor  response 
and  for  the  imaging  approach  was  obtained  through  testing  over  known  objects.  These  tests 
included  measurements  over  aluminum  spheres  and  known  UXO. 

Initial  testing  was  performed  on  an  assortment  of  aluminum  spheres  ranging  from  3.8  to  20.3  cm 
(1.5  to  8  in.)  in  diameter  in  addition  to  representative  munitions  and  clutter.  Since  an  existing 
drive  winding  was  used  in  this  program,  which  has  a  depth  limitation  of  roughly  0.45  m  (18 
in.),  including  the  standoff  distance  to  the  ground  surface,  the  focus  has  been  on  objects 
located  within  approximately  0.4  m  (16  in.)  of  the  surface.  Since  only  the  smaller  UXO  were 
within  this  depth  at  the  Aberdeen  Proving  Grounds  Test  Site,  only  the  smaller  UXO  were 
ordered.  A  new  drive  is  to  be  designed  in  the  follow-on  effort  to  extend  depth  of  sensitivity  to 
greater  than  0.6  m  (24  in.)  with  a  goal  of  a  0.9  m  (36  in.)  depth  of  sensitivity.  (Note  that  the 
object  or  sphere  size  has  not  yet  been  determined  for  this  requirement.) 

Figure  12  shows  representative  scans  of  objects  over  the  older  version  of  the  array  (with  8 
sensing  elements).  These  scans  simulated  a  scan  rate  of  approximately  1.2  m/sec  (4  ft/sec)  and  a 
distance  to  the  array  of  15  cm  (6  in.).  The  bomblet  appears  as  a  discrete  object  while  the 
dimensions  of  the  cylinder  are  large  compared  to  the  width  of  the  array.  Even  in  this  simple 
case,  the  image  readily  shows  the  shape  and  differences  for  the  two  objects  and  the  orientation 
variation  for  the  cylinder.  The  non-uniformity  of  the  cylinder  image  is  attributed  to  the 
variations  in  the  magnetic  field  intensity  across  the  array.  This  can  be  corrected  with  model- 
based  methods  that  correctly  account  for  the  field  variations  at  the  sides  of  the  drive  winding,  by 
widening  the  drive  winding,  or  by  improving  the  drive  winding  designs  to  reduce  this  field 
variation. 

Figure  13  shows  the  data  and  measurement  grids  from  one  of  the  channels  beneath  the  object. 
This  grid  was  generated  for  the  bomblet  using  the  revised  sphere  model  and  illustrates  how  the 
bomblet  properties  (size  and  depth)  can  be  estimated  from  the  grid.  The  cylinder  data  does  not 
even  fall  on  the  grid.  This  illustrates  the  basic  capability  for  object  identification/discrimination 
and  clutter  suppression  but  more  robust  procedures  still  need  to  be  developed.  For  example, 
different  measurement  grids  need  to  be  used  for  the  different  sensing  elements  to  account  for  the 
variations  in  the  magnetic  field  intensity  across  the  array. 
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Figure  12.  Representative  images  of  data  taken  with  the  old  array. 


Figure  13.  Representative  measurement  grid,  data,  and  depth/size  estimates  for  a  bomblet  and  cylinder  scanned 
across  the  “old”  sense  coil. 


After  fabrication  of  the  higher  resolution  array  of  sense  elements  and  integration  with  the 
impedance  measurement  instrumentation,  tests  were  performed  on  single  objects  at  various 
locations  beneath  the  sensor  array.  Figure  14  shows  the  basic  configuration  for  the  array  along 
with  the  numbering  for  the  sense  elements  and  the  locations  of  the  aluminum  spheres.  Tests 
were  also  performed  with  0.9  m  (36  in.)  long  pipes  and  several  UXO. 

Figure  15  shows  the  response  of  several  sense  elements  as  different  size  aluminum  spheres  are 
lowered  beneath  sense  element  8.  The  model  was  used  to  calculate  an  appropriate  measurement 
grid  for  the  sphere  location  and  properties  (e  g.,  size  and  electrical  conductivity)  and  for  the 
sense  element  position  and  orientation  within  the  drive  winding.  When  the  sense  element  is 
directly  over  the  object,  in  this  case  element  8,  the  grid  indicates  reasonable  sensitivity  to  the  size 
and  depth  of  the  object,  as  indicated  by  the  relatively  large  spacing  between  the  grid  lines.  For 
each  sphere,  the  model  and  data  show  a  nonlinear  response  and  the  data  follows  a  line  of 
constant  size  as  the  depth  is  varied.  The  response  of  offset  sense  elements,  in  this  case  elements 
4  and  12,  show  similar  behavior  to  sense  element  8,  but  with  a  reduced  sensitivity  to  the  object 
properties.  This  indicates  that  the  adjacent  sense  elements  provide  complementary  information 
about  the  object  properties.  The  data  does  not  fall  on  the  grids  for  these  preliminary 


Page  27  of  66 


measurements  because  the  sphere  location  beneath  the  sense  elements  was  only  controlled  by 
hand,  the  thickness  of  the  sphere  walls  is  not  accounted  for  in  the  model  and  the  model  does  not 
currently  account  for  the  field  variations  over  the  footprint  of  the  sense  elements.  However,  even 
so,  the  data  is  relatively  close  to  the  grid  and  follows  predicted  behavior  as  lift-off  is  varied. 

In  this  case,  the  measurement  grid  is  plotted  as  the  imaginary  part  of  the  signal  against  the  real 
part  of  the  signal.  This  is  more  convenient  (in  some  cases)  than  the  magnitude  and  phase  for 
data  plotting,  particularly  when  the  measurement  signal  is  small  and  comparable  to  the 
measurement  noise.  This  is  because  the  drive  winding  enhancement  shifts  the  air  point  to  the 
origin  and  measurement  noise,  particularly  near  the  air  point,  can  give  rise  to  any  phase  angle 
between  -180  and  +180  degrees.  The  result  is  a  spread  of  data  points  over  all  phase  angles,  for 
small  magnitudes,  instead  of  a  clustering  of  data  points  near  the  origin.  Plots  of  the  magnitude 
itself  are  still  convenient,  though,  since  it  combines  the  real  and  imaginary  parts  of  the  signal  into 
a  single  value  that  can  be  displayed. 

Figure  16  shows  the  estimated  properties  obtained  from  sense  element  8  as  aluminum  spheres  of 
different  size  are  lower  beneath  it.  The  depth  estimates  are  generally  within  approximately  2  cm 
to  a  depth  of  approximately  20  cm,  which  is  where  the  noise  in  the  measurements  causes  some  of 
the  data  to  fall  off  of  the  measurement  grid  and  the  standard  deviation  of  the  measurement, 
indicated  by  the  vertical  error  bars,  increases  substantially.  The  associated  estimates  of  the 
sphere  sizes  are  in  reasonable  agreement  with  the  known  sphere  sizes  for  shallow  depths.  These 
estimates  become  less  reliable  at  depths  of  around  20  cm  because  of  the  decay  of  the  imposed 
magnetic  field  intensity  with  depth.  Use  of  a  larger  drive  winding  wavelength  would  reduce  this 
decay,  extending  the  field  depth  of  penetration.  Also,  better  size  estimates  at  deeper  depths 
possible  with  the  larger  objects.  Figure  17  shows  the  sensor  response  and  measurement  grids  for 
spheres  beneath  sense  elements  that  are  closer  to  the  edges  of  the  array.  Similar  to  Figure  15, 
there  is  good  sensitivity  to  both  the  sphere  size  and  depth. 

Figure  18  shows  the  effect  of  an  offset  of  the  sphere  on  the  measurement  response.  This  would 
typically  result  from  a  scan  of  an  array  over  an  object  or  when  sense  elements  are  placed  at 
various  locations  within  the  footprint  of  the  drive  winding  and  exposed  to  different  magnetic 
field  intensities.  In  this  case,  the  aluminum  spheres  are  offset  by  15.2  cm  (6  in.)  on  either  side  of 
the  linear  row  of  sense  elements  ( y  direction),  but  aligned  with  sense  element  8  in  the  x  direction. 
The  measurement  grids,  which  account  for  the  offset  in  object  and  sense  element  positions, 
indicate  that  the  response  should  be  symmetric,  since  the  sense  elements  are  placed  at  the  center 
of  the  drive  winding  construct.  However,  the  data  is  not  symmetric,  indicating  differences  in 
either  the  actual  positions  of  the  spheres  (e.g.,  misalignment  with  element  8)  or  an  unmodeled 
asymmetry  in  the  drive  winding. 

This  shows  that  the  model  captures  the  basic  response  of  the  sensors  to  spheres  under  the 
footprint  of  the  array  and  with  various  positions  of  the  sphere  relative  to  the  sense  elements. 

Each  of  the  sense  elements  provide  information  about  the  object  properties;  the  capability  for 
using  this  information  is  related  to  the  sensitivity,  which  can  be  determined  from  the  model,  and 
the  measurement  noise,  which  is  determined  empirically.  In  follow-on  efforts  more  systematic 
methods  will  be  applied  for  using  information  from  multiple  sensing  elements  to  estimate 
multiple  properties  of  objects  such  as  size,  depth,  shape  and  orientation. 
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Figure  14.  Schematic  for  single  object  measurements  with  only  depth  varied. 
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Figure  15.  Response  of  a  single  sense  element  (number  4,  8,  or  12)  as  the  depth  of  several  aluminum  spheres, 
placed  beneath  sense  element  8,  is  varied. 
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Figure  16.  Estimated  depth  and  size  from  element  8  for  aluminum  spheres  placed  beneath  the  element 
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Figure  17.  Response  of  an  element  as  the  depth  of  several  aluminum  spheres  beneath  the  element  is  varied. 


Figure  18.  Response  of  sense  element  8  as  the  depth  of  several  aluminum  spheres  is  varied.  The  spheres  are  aligned 
with  element  8  but  offset  by  6  in.  across  the  drive  windings. 
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Measurements  of  objects  beneath  a  sense  element  as  the  object  depth  is  varied  also  indicate  that 
the  real  and  imaginary  parts  of  the  signal  can  be  used  to  discriminate  between  objects.  Figure  19 
shows  the  response  as  the  depth  of  several  UXO  is  varied  from  approximately  3.8  cm  (1.5  in.) 
from  the  array  to  the  top  of  the  object,  to  greater  than  30.5  cm  (12  in.).  Except  for  the  spherical 
BLU-26  the  response  was  measured  for  several  orientations  of  the  major  axis  (longest 
dimension)  either  parallel  or  perpendicular  to  the  linear  array  of  sense  elements,  or  vertical  with 
the  UXO  pointed  down  (tail  up)  or  up  (tail  down).  In  each  case,  the  response  converges  to  the 
air  point  when  the  object  is  sufficiently  far  from  the  array  to  be  detected.  For  the  M385,  all  three 
orientations  provide  essentially  the  same  response,  with  a  large  real  component  and  a  small 
imaginary  component.  For  the  BLU-26,  the  real  and  imaginary  components  are  comparable. 

For  the  higher  aspect  ratio  M86  and  the  MK2,  the  parallel  and  perpendicular  orientations  are 
similar  but  the  vertical  orientation  provides  a  noticeable  shift  in  the  response.  Note  that  the 
vertical  M86  was  pointed  up  in  an  attempt  to  better  align  the  UXO  with  the  sense  element.  The 
response  shift  is  larger  for  the  larger  M86  UXO.  Some  of  the  reduction  in  signal  associated  with 
the  vertical  orientation  is  due  to  the  center  of  the  object  being  effectively  deeper  for  the  vertical 
orientation.  The  similarities  in  the  responses  for  the  M86  and  MK2,  at  least  in  the  parallel  and 
perpendicular  orientations,  along  the  depth  variation  line  indicates  potential  non-uniqueness 
issues  for  a  single  sense  element  and  a  simple  sphere  model  response;  the  objects  give  a  similar 
response  at  different  depths,  even  though  they  are  different  sizes.  Other  information,  for 
example  from  other  sense  elements  (e  g.,  images),  other  sense  element  orientations,  or  multiple 
frequencies,  would  be  needed  to  properly  discriminate  between  these  objects.  This  further 
supports  the  value  of  a  two-dimensional  array  of  sensing  elements. 

Additional  tests  were  performed  by  scanning  the  array  over  known  objects  to  demonstrate  the 
imaging  capability  of  the  array.  In  the  laboratory,  this  was  performed  using  the  test 
configuration  of  Figure  6,  which  has  the  sensor  array  mounted  within  a  cart  structure.  A  rotary 
position  encoder  provides  the  linear  position  of  the  array  along  the  scan  path.  Objects  were  then 
placed  at  a  fixed  location  and  orientation  beneath  the  array  while  a  scan  was  performed. 


Figure  19.  Response  of  sense  element  8  as  the  depth  of  several  objects  is  varied. 
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As  a  preliminary  baseline  measurement,  a  scan  was  taken  without  any  objects  beneath  the  array. 
The  results  are  shown  in  Figure  20  and  indicate  modest  variation  of  the  sensor  response  between 
the  elements.  Ideally  the  response  is  zero  in  air  and  any  variations  in  the  response  can  be 
attributed  to  the  instrument  noise.  Vertical  stripes  (common  to  all  sense  elements)  corresponds 
to  effects  from  the  cart  motion  while  the  gradual  variation  in  response  along  the  length  of  the 
scan  corresponds  to  drift  in  the  response.  The  array  was  calibrated  just  prior  to  performing  this 
scan.  Figure  21  provides  a  comparison  plot  between  the  sense  elements.  Note  that  the  units  for 
the  plots  are  different  than  the  units  for  the  images,  which  were  obtained  directly  from  the 
GridStation  software,  because  the  signal  from  the  coils  was  unconventionally  calibrated  to  be  in 
terms  of  the  field  intensity  instead  of  a  voltage. 


Figure  20.  Array  magnitude  response  in  air,  approximately  1.2  m  (48  in.)  above  a  cement  floor.  The  horizontal  axis 
gives  the  scan  distance  in  inches.  The  vertical  axis  gives  the  distance  across  the  sense  element  array  in  inches. 


Figure  21.  Plots  of  the  average  response  of  each  element  of  the  array  for  the  scan  in  air  of  Figure  20.  The  error  bars 
provide  the  standard  deviation  of  the  measurement. 
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Representative  scan  images  over  a  12.7  cm  (5  in.)  diameter  aluminum  sphere  are  shown  in 
Figure  22  for  a  5  cm  (2  in.)  depth.  This  depth  is  given  to  the  top  of  the  sphere.  The  figure  shows 
images  of  the  raw  data  in  the  form  of  the  magnitude,  real  part,  and  imaginary  part  of  the  signal. 

In  this  case  the  real  part  of  the  signal  is  substantially  larger  than  the  imaginary  part  so  that  the 
magnitude  tends  to  reflect  the  real  part  of  the  signal.  The  sphere  is  readily  discerned  from  any 
background  variations. 
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Figure  22.  Array  response  over  a  12.7  cm  (5  in.)  diameter  aluminum  sphere  at  a  depth  of  5  cm  (2  in.),  (a) 
Magnitude,  (b)  Real  part,  (c)  Imaginary  part  The  horizontal  axis  gives  the  scan  distance  in  inches.  The  vertical 
axis  gives  the  distance  across  the  sense  element  array  in  inches. 


These  measurements  also  demonstrate  the  diffusion  or  spreading  of  the  field  response  as  objects 
get  farther  away  from  the  sensor  array.  The  deeper  sphere  has  a  larger  effective  spatial  image 
than  the  closer  objects.  This  is  illustrated  clearly  in  Figure  23  for  images  taken  for  several  depths 
of  the  12.7  cm  (5  in.)  diameter  sphere.  The  contrast  is  also  evident  in  Figure  25,  where  the 
smaller  but  deeper  sphere  provides  a  larger  spatial  image  than  the  larger  but  shallower  sphere. 
Thus,  the  image  by  itself  cannot  be  used  to  determine  the  object  size  and  other  information,  such 
as  the  magnitude  scaling,  is  necessary  to  correctly  differentiate  the  object  properties.  If  the  depth 
is  known,  then  the  image  size  directly  reflects  the  object  size,  as  shown  in  Figure  26  for  several 
different  sized  spheres  at  a  fixed  depth,  which  is  taken  as  the  distance  to  the  top  of  the  spheres. 
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The  measurement  grid  can  be  used  in  conjunction  with  the  image  scans  to  determine  effective 
properties  for  hidden  objects.  This  is  illustrated  in  Figure  24  for  the  12.7  cm  (5  in.)  diameter 
aluminum  sphere.  The  data  points  were  taken  from  the  peaks  in  the  images  in  the  real  and 
imaginary  components  of  the  impedance.  Note  that  Figure  23  only  shows  the  magnitude 
response.  This  data  compares  favorably  with  the  measurements  taken  as  the  aluminum  sphere 
was  lowered  by  hand  beneath  sense  element  8.  Differences  in  the  measurements  can  be 
attributed  to  hand  held  measurements  being  slightly  off  the  center  line  of  the  sense  element. 


Figure  23.  Array  magnitude  response  over  a  12.7  cm  (5  in.)  diameter  aluminum  sphere  at  a  depth  of  (a)  5  cm  (2  ia), 
(b)  10  cm  (4  in.),  (c)  15  cm  (6  in.),  and  (d)  20  cm  (8  in.). 
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Figure  24.  Response  of  the  sensor  array  to  the  12.7  cm  (5  in)  diameter  aluminum  sphere.  The  grid  is  for  sense 
element  8.  The  stationary  cart  data,  where  the  sphere  was  held  beneath  element  8  and  lowered  by  hand  is  plotted 
along  with  impedance  values  taken  from  the  peak  signal  responses  in  the  image  scans. 
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Figure  25.  Array  magnitude  response  over  aluminum  spheres  of  (a)  12.7  cm  (5  in.)  diameter  and  a  depth  of  5  cm  (2 
in.)  and  (b)  6.4  cm  (2.5  in.)  diameter  and  a  depth  of  20  cm  (8  in). 
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Figure  26.  Array  magnitude  response  over  aluminum  spheres  at  a  depth  of  10  cm  (4  in.)  for  diameters  of  (a)  6.4  cm 
(2.5  in.),  (b)  7.6  cm  (3  in.),  (c)  10.2  cm  (4  in.),  and  (d)  12.7  cm  (5  in.). 
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Scan  images  from  the  array  can  also  reflect  the  orientation  of  hidden  objects.  To  demonstrate 
this  effect,  two  0.9  m  (36  in.)  long  pipes  were  placed  at  various  orientations  and  scanned  with  the 
array.  A  representative  set  of  measurement  images  is  shown  in  Figure  27  for  a  63.8  cm  (1.5  in.) 
diameter  pipe  oriented  parallel  to  the  drive  winding  segments  and  perpendicular  to  the  scan 
direction.  The  impedance  data  (magnitude,  real  part,  and  imaginary  part)  shows  an 
approximately  equal  response  for  each  of  the  array  elements,  except  for  the  outermost  elements. 
This  indicates  that  each  sense  element  sees  essentially  the  same  magnetic  field;  the  outermost 
sense  elements  are  most  likely  influenced  by  the  return  legs  of  the  drive  windings.  Extending 
the  drive  winding  segments  further  past  the  sense  element  array  and  modifying  the  drive 
winding  design  to  reduce  the  effect  of  these  return  segments  should  improve  the  quality  of 
the  images. 
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Figure  27.  Array  response  over  a  91  cm  (36  in.)  long  3.8  cm  (1.5  in.)  diameter  pipe  oriented 
parallel  to  the  drive  windings  and  at  a  depth  of  10  cm  (4  in.),  (a)  Magnitude,  (b)  Real  part,  (c) 
Imaginary  part.  The  horizontal  axis  gives  the  scan  distance  in  inches.  The  vertical  axis  gives 
the  distance  across  the  sense  element  array  in  inches. 
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Similar  results  are  obtained  with  other  orientations  and  other  depths.  For  example,  Figure  28 
shows  the  magnitude  response  for  three  different  orientations  of  the  pipe:  parallel,  perpendicular, 
and  at  45  degrees  to  the  sense  element  array.  The  parallel  and  45  degree  orientation  images 
show  a  reasonably  uniform  response  except  for  the  end  effects  at  the  edges  of  the  array.  Figure 
29  shows  the  images  for  a  deeper  pipe.  Note  again  that  the  spatial  width  of  the  images  increases 
with  object  depth,  corresponding  to  the  greater  diffusion  of  the  magnetic  fields  from  the  deeper 
object.  Figure  30  and  Figure  3 1  show  the  corresponding  data  for  the  larger  diameter  pipe  at  two 
depths.  The  width  of  the  detected  object  is  larger,  as  expected,  since  the  pipe  is  larger  in  this 
case  while  the  depths  were  kept  the  same. 


■  I  I  I  . i  .  '  '  "i"7—  i  I 

0.1  0.2  0.3  0.4  05  as  07  08  09  1.0  1.1  1.2  1.3 

(uH) 

(C) 

Figure  28.  Array  magnitude  response  over  a  91  cm  (36  in.)  long  3.8  cm  (1.5  in.)  diameter  pipe  at  a  depth  of  10  cm 
(4  in.)  oriented  (a)  parallel  to  the  sense  element  array,  (b)  at  45  degrees  to  the  sense  element  array,  and  (c) 
perpendicular  to  the  sense  element  array. 
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Figure  29.  Array  magnitude  response  over  a  91  cm  (36  in.)  long  3.8  cm  (1.5  in.)  diameter  pipe  at  a  depth  of  20  cm 
(8  in.)  oriented  (a)  parallel  to  the  sense  element  array,  (b)  at  45  degrees  to  the  sense  element  array,  and  (c) 
perpendicular  to  the  sense  element  array. 
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Figure  30.  Array  magnitude  response  over  a  91  cm  (36  in.)  long  5. 1  cm  (2  in.)  diameter  pipe  at  a  depth  of  10  cm  (4 
in.)  oriented  (a)  parallel  to  the  sense  element  array,  (b)  at  45  degrees  to  the  sense  element  array,  and  (c) 
perpendicular  to  the  sense  element  array. 
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Figure  3 1.  Array  magnitude  response  over  a  91  cm  (36  in.)  long  5. 1  cm  (2.0  in.)  diameter  pipe  at  a  depth  of  20  cm 
(8  in.)  oriented  (a)  parallel  to  the  sense  element  array,  (b)  at  45  degrees  to  the  sense  element  array,  and  (c) 
perpendicular  to  the  sense  element  array. 


Figure  32  through  Figure  35  show  some  scan  images  over  several  UXO  having  different 
orientations  at  a  fixed  depth.  In  this  images,  the  depth  is  the  distance  from  the  sensor  array  to  the 
top  of  the  UXO.  Figure  32  shows  the  basic  response  for  a  bomblet  which  does  not  change 
substantially  with  orientaions.  The  other  UXO  (MK2  for  Figure  33,  M86  for  Figure  34,  and 
MK1 18  Rockeye  for  Figure  35)  have  a  longer  axis  that  appears  in  the  45°  angle  images  and  is 
most  evident  when  the  UXO  is  perpendicular  to  the  array  of  sense  elements.  The  parallel 
orientation  does  not  show  an  elongation  of  the  image  because  the  end  effects  at  the  edges  of  the 
array  cause  an  effective  narrowing  of  the  object  image.  This  narrowing  effect  is  illustrated  in  the 
parallel  orientation  images  over  long  pipes  of  Figure  27  through  Figure  31.  Extending  the 
lengths  of  the  drive  winding  segments  and  using  other  drive  winding  configurations  that  reduce 
the  inhomogeneity  across  the  sense  element  array  should  provide  images  that  show  the  object 
orientation  in  the  parallel  direction.  Alternative,  compensation  techniques  can  be  applied  to  the 
images  themselves,  as  illustrated  in  Figure  36,  to  correct  for  unmodeled  inhomgeneneities. 
Another  difference  is  the  spatial  resolution  of  the  image.  In  the  scan  direction,  the  increment 
between  data  points  is  typically  less  than  1  mm  (0.039  in.).  In  contrast,  in  the  array  direction,  the 
distance  between  elements  is  25.4  mm  (1.0  in.)  and  the  image  resolution  is  lower.  Thus,  use  of 
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Figure  33.  Array  magnitude  response  over  an  MK2  for  several  orientations  with  respect  to  the  sense  element  array 
with  the  distance  from  the  array  to  the  top  of  the  UXO  10.2  cm  (4  in).  The  axes  give  the  scan  dimensions  in  inches. 
The  dashes  indicate  a  circular  response 
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Figure  35.  Array  magnitude  response  over  an  MK1 18  Rockeye  for  several  orientations  with  respect  to  the  sense 
element  array  with  the  distance  from  the  array  to  the  top  of  the  UXO  10.2  cm  (4  in.).  The  axes  give  the  scan 
dimensions  in  inches.  The  dashes  indicate  a  circular  response. 
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Figure  36.  Compensated  array  magnitude  response  over  an  MK2  for  several  orientations  which  correspond  to 
Figure  33. 


The  magnitude  images  of  Figure  32  through  Figure  35  have  embedded  within  them  the 
impedance  images  and  data  that  can  provide  information  about  the  object  characteristics.  The 
peak  values  of  the  real  and  imaginary  portion  of  the  response  were  taken  from  the  scan  images 
and  plotted  in  Figure  37.  The  lines  indicate  the  response  of  the  sense  element  if  the  sense 
element  was  centered  and  raised  over  the  object  from  Figure  19.  The  different  UXO  have 
different  base  impedance  values,  consistent  with  the  measurements  taken  as  the  UXO  are 
dropped  below  a  specific  sense  element,  as  in  Figure  19.  The  only  significant  difference  in 
responses  occurs  for  the  vertical  M86,  which  was  pointing  upward  in  Figure  19  and  downward  in 
Figure  37.  Note  that  scan  images  also  tend  to  have  slightly  smaller  values  for  the  peak  response 
because  the  image  is  providing  a  spatial  average  of  values  in  two  directions  (x  andy)  while  the 
single  element  measurement  do  not  have  any  spatial  averaging  of  the  measurement  signal.  This 
confirms  that  the  scan  images  can  provide  some  impedance  information  for  providing  object 
discrimination. 

These  values  for  the  measured  response  were  then  processed  using  sphere  model  grids  to 
determine  the  effective  size  and  depth  of  the  object.  The  same  grid  was  used  for  each  orientation 
of  an  object.  Different  grids,  based  on  nominal  electrical  properties,  were  used  for  each  object, 
except  for  the  MK2  and  the  M86  which  used  the  same  grid.  These  results,  given  in  Table  1,  are 
in  reasonable  agreement  with  the  nominal  object  sizes.  Note  that  the  UXO  are  generally  not 
spherical  so  that  the  sphere  model  is  not  truly  appropriate;  however,  it  does  provide  a  reasonable 
estimate  of  the  object  properties.  Furthermore,  this  is  a  simple  use  of  a  model  for  estimating 
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properties.  Only  a  single  value  from  the  scan  image  is  used  to  estimate  the  object  properties. 

Better  estimates  are  expected  when  the  nonspherical  shape  is  accounted  for  and  the 
responses  from  multiple  sense  elements  are  used.  Thus,  the  use  of  ellipsoid  models, 
multiple  sensing  elements  to  discriminate  object  features  (shape,  size,  orientation)  should 
be  investigated  in  the  follow-on  effort. 


Figure  37.  Impedance  response  for  several  UXO  taken  from  scan  images  with  the  distance  from  the  array  to  the  top 
of  the  UXO  10.2  cm  (4  in.). 


for  the  UXO  measurements  of  Figure  37. 
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Measurements  have  also  been  performed  over  multiple  objects  in  order  to  show  the  benefit  of 
high  resolution  imaging.  Figure  38  shows  several  response  images  for  the  bomblets  aligned 
parallel  to  the  sense  element  array.  Figure  39  shows  several  response  images  for  the  bomblets 
aligned  perpendicular  to  the  sense  element  array.  At  shallow  depths  and  with  large  separation 
distances,  the  response  of  each  bomblet  is  distinct.  For  deeper  depths  and  smaller  separation 
distances,  the  interactions  between  the  objects  become  more  apparent  and  the  diffusion  of  the 
magnetic  field  response  leads  to  an  overlap  of  the  effective  response  and  an  apparent  “single” 
object  response.  This  response  has  an  effective  orientation  associated  with  alignment  of  the 
objects  and  also  local  maxima  that  are  associated  with  the  center  of  the  objects.  For  very  deep 
objects  and  for  very  small  separations,  the  effective  orientation  of  the  response  becomes  less 
apparent  and  symmetric  responses  result.  The  use  of  a  high  resolution  imaging  array  allows  the 
spatial  variations  of  the  response  from  the  objects  to  be  captured,  including  orientation  and  local 
maxima  and  minima  in  the  response  that  may  permit  separation  of  object  responses.  Combining 
information  from  high  resolution  images  with  property  estimates  from  measurement  grids 
at  multiple  sensing  elements  should  be  addressed  in  the  follow-on  effort. 


(uH) 

Figure  38.  Array  response  over  two  BLU-26  bomblets  aligned  parallel  to  the  sense  element  array.  The  axes  give 
the  image  dimensions  in  inches. 
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Figure  39.  Array  response  over  two  BLU-26  bomblets  10.2  cm  (4  in.)  below  the  array.  The 
inside  spacing  between  the  bomblets  is  indicated.  The  axes  give  the  image  dimensions  in 
inches. 
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Field  Tests 

Preliminary  tests  were  performed  on  the  calibration  grid  of  the  Aberdeen  Proving  Grounds 
Standardized  Test  Site  in  December,  2002.  The  goal  of  these  measurements  was  to  obtain  initial 
field  experience  at  a  test  site,  which  would  provide  input  for  enhancement  criteria  and 
recommendations.  The  focus  of  these  tests  was  on  the  objects  located  nearer  the  surface,  with  a 
depth  less  than  0.4  m,  which  was  generally  the  smaller  objects. 

A  schematic  diagram  of  the  objects  and  their  locations  on  the  test  grid  is  shown  in  Figure  41 . 

The  test  grid  is  divided  into  1  m  sections  with  letters  denoting  the  grid  column  and  numbers 
denoting  the  grid  lane.  An  object  is  placed  near  the  center  of  various  grid  cells  with  depths  and 
orientations  varying.  Eight  pound  shot  put  markers  are  placed  around  the  edge  of  the  calibration 
grid  and  at  several  locations  within  the  grid.  The  distance  between  objects  is  generally  2  m  or 
greater. 

40  n - 1 - 1 - i - 1 - 1 - n 


30 

"E 


E 

o 

Z 

10 


0 3  -e* <* — m — *, — & - — ---■i-t 

0  5  10  15  20  25  30 

Nominal  North  (m) 

Figure  40.  Schematic  diagram  of  the  calibration  grid. 
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Measurements  were  performed  at  various  locations  in  the  grid  using  the  instrument  configuration 
shown  in  Figure  43 .  This  was  essentially  the  same  as  the  configuration  used  in  the  laboratory. 
Power  was  supplied  by  a  portable  generator  located  at  a  distance  from  the  equipment,  near  the 
vehicles.  The  original  intent  for  the  measurements  was  to  position  the  table  with  the  support 
instrumentation  (impedance  instrument  and  laptop)  near  several  grid  cells.  Then,  measurements 
would  be  performed  by  scanning  the  cart  over  these  cells  without  moving  the  table.  However,  it 
was  decided  in  the  field  that  it  would  be  more  convenient  for  covering  larger  areas  to  simply 
move  the  table  along  with  the  cart  so  that  longer  measurement  scans  could  be  obtained.  As  a 
result,  for  measurement  scans  covering  distances  less  than  approximately  2  m  the  table  was  not 
moved,  but  for  longer  scans,  the  table  was  moved  along  side  of  the  cart.  Scan  distances  were 
monitored  with  a  rotary  position  encoder  on  the  rear  wheel  axel  of  the  cart.  The  encoder  was 
calibrated  with  a  tape  measured  and  periodically  checked  during  the  trial. 
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Figure  41.  Photograph  of  the  cart  with  sensor  array  and  support  instrumentation  used  during  the  field  test. 


The  site  was  relatively  wet  during  these  tests.  The  week  prior  to  these  field  tests,  they  had 
received  several  inches  of  snow,  which  was  followed  by  several  days  of  rain.  The  ground  had  a 
day  or  two  to  dry  out  slightly  before  our  arrival.  During  the  testing  the  mornings  were  cold 
enough  that  the  ground  was  basically  frozen,  but  in  the  afternoons  it  had  warmed  enough  so  that 
some  areas  had  free  standing  water  and  the  ground  was  a  relatively  soft  mud.  Fortunately,  most 
of  the  shallower  objects  of  interest  were  located  in  lanes  1  through  8  which  were  on  relatively 
dryer  and  flat  ground.  Figure  41  shows  that  the  ground  was  relatively  flat  in  some  areas.  Some 
of  the  other  objects,  particularly  those  in  the  lower  left  comer  of  the  grid  were  under  water  or  in 
relatively  rough  terrain,  as  shown  in  Figure  42.  When  possible,  attempts  were  made  to  scan  over 
the  objects  in  these  areas. 


Figure  42.  Photographs  of  some  of  the  areas  of  rough  terrain,  some  of  which  were  also  under  water. 


Measurements  were  generally  performed  in  a  scanning  mode,  where  the  cart  was  pushed  along 
the  test  lane.  Figure  43  provides  a  schematic  for  the  scan  of  the  cart  along  the  center  of  a  path. 
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The  sense  elements  only  span  a  portion  of  the  cell  so  complete  coverage  of  the  cell  was  obtained 
with  multiple  scan  passes.  Typical  offsets  of  the  array  were  approximately  0.36  m  (14  in.)  on 
either  side  of  the  center  scan  path  width  from  the  center  so  that  there  were  several  elements  of 
overlap  with  the  “center”  scan  path.  The  positions  of  the  grid  cell  markers  were  confirmed 
visually  with  the  row  and  column  flags  and  the  other  cell  markers.  They  were  also  confirmed 
with  a  tape  measure  to  be  in  approximately  the  correct  location. 


© 


Grid  Cell 
Markers 


Figure  43.  Schematic  diagram  for  the  cart  and  sensor  array  orientation  for  a  scan  along  the  center  of  a  grid  cell. 


For  these  scans  over  the  ground  surface,  the  typical  standoff  distance  or  lift-off  between  the 
sensor  array  and  the  ground  surface  was  about  5  cm  (2  in.)  but  was  increased  or  decreased  as 
indicated  for  individual  scans.  The  nominal  distance  changed  as  the  array  was  scanned  over  the 
surface  since  the  ground  was  rough,  with  some  deep  ruts  that  led  in  some  case  to  ~  10  cm  (4  in.) 
air  gaps  and  also  caused  some  tilt  in  the  cart  holding  the  array.  There  was  also  some  settling  in 
the  areas  where  objects  were  buried. 

A  representative  “vertical”  scan  down  a  column  is  shown  in  Figure  44.  As  expected  from  the 
laboratory  measurements,  the  different  objects  had  distinguishable  signature  responses. 
The  M55  appears  in  the  imaginary  part  image,  but  not  in  the  real  part  image.  In  contrast,  the 
M3 85  appears  predominantly  in  the  real  part  image,  with  a  relatively  weak  response  in  the 
imaginary  part  image.  The  BLU-26  provides  a  strong  response  in  both  the  real  and  imaginary 
part  images,  while  the  deep  M86  is  not  apparent.  When  only  the  magnitude  of  the  response  is 
plotted,  which  in  this  case  is  dominated  by  the  real  part  of  the  response,  the  indication  for  the 
M55  is  not  obvious.  However,  a  logarithmic  color  scale  accentuates  the  low  signals  and 
compresses  the  high  signals,  which  allows  the  M55  to  be  revealed  in  the  magnitude  image,  as 
shown  in  Figure  45.  A  drawback  is  that  more  of  the  noise  is  also  visible  in  the  image.  For 
convenience  in  plotting  the  large  number  of  scans  over  the  test  grids,  much  of  the  scan  data  will 
be  presented  as  magnitude  images  using  the  logarithmic  color  scale.  The  repeatability  of  the 
scans  is  demonstrated  in  Figure  46. 
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Figure  44.  Vertical  scan  down  the  center  of  cells  A3-A6.  Tie  horizontal  axis  gives  the  scan  distance  in  inches.  The 
vertical  axis  gives  the  distance  along  the  sense  element  array  and  spans  0.41  m  (16  in.).  The  object  type,  depth, 
azimuth  angle,  dip  angle  and  approximate  position  are:  M55  (0. 1  m,  0°,  45°)  at  -28  in.,  M385  (0. 1  m,  0°,  45°)  at 
-106  in.,  BLU-26  (0. 1  m,  0°,  0°)  at  -184  in.,  and  M86  (0.4  m,  0°,  45°)  at  -262  in.  The  top  image  gives  the 
magnitude,  the  middle  gives  the  real  part,  and  the  bottom  gives  the  imaginary  part 


Figure  45.  A  plot  of  the  magnitude  response  from  Figure  44  using  a  logarithmic  color  scale. 


Figure  46.  Repeatability  of  magnitude  images  for  vertical  scans  down  the  center  of  cells  A3-A6.  The  horizontal 
axis  gives  the  scan  distance  in  inches.  The  vertical  axis  gives  the  distance  along  the  sense  element  array  and  spans 
0.4 1  m  (16  in).  The  object  type,  depth,  azimuth  angle,  dip  angle  and  approximate  position  are:  M55  (0. 1  m,  0°,  45°) 
at  -28  in.,  M385  (0.1  m,  0°,  45°)  at  -106  in.,  BLU-26  (0.1  m,  0°,  0°)  at  -184  in.,  and  M86  (0.4  m,  0°,  45°)  at -262  in. 
A  calibration  was  only  performed  prior  to  the  first  scan. 
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Similar  measurement  scans  were  performed  with  the  array  offset  from  the  center  of  the  objects, 
as  shown  in  Figure  47.  In  this  case  the  array  was  moved  in  approximately  0. 1 5  m  (6  in.) 
increments  from  the  left  side  of  the  grid  ceils  to  the  right  side  of  the  grid  cells.  Again,  the  M3  85 
and  BLU-26  are  distinct,  particularly  when  they  are  near  the  center  of  the  array.  When  they  are 
underneath  the  outermost  (left  or  right  side)  sense  elements,  edge  effects  due  to  the  proximity  of 
the  sides  or  return  legs  of  the  drive  winding  distort  the  object  images. 


0.3  m  left  shift  0.15  m  left  shift  middle  scan  0.15  mright  shift  0.3  m right  shift 

Figure  47.  Effect  of  array  offset  on  image  response  for  vertical  scans  down  cells  A3-A6.  The  vertical  axis  gives  the 
scan  distance  in  inches.  The  horizontal  axis  gives  the  distance  along  the  sense  element  array  and  spans  0.41  m  (16 
in.).  The  object  type,  depth,  azimuth  angle,  dip  angle  ami  approximate  position  axe:  M55  (0. 1  m,  0°,  45°)  at  ~28  in., 
M385  (0. 1  m,  0°,  45°)  at  -106  in.,  BLU-26  (0. 1  m,  0°,  0°)  at  -184  in.,  and  M86  (0.4  m,  0°,  45°)  at  -262  in. 
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Measurements  were  also  performed  at  several  different  lift-offs  or  stand-off  distances  from  the 
ground.  Images  of  the  real  and  imaginary  parts  of  the  response,  shown  in  Figure  48,  show  the 
decreasing  sensitivity  to  the  presence  of  the  objects  as  the  array  is  raised  higher  above  the 
ground.  In  this  case,  the  objects  were  nominally  10  cm  beneath  the  surface  except  for  the  M86, 
which  was  buried  to  a  depth  of  40  cm.  Note  that  the  nominal  stand-off  distance  was  5. 1  cm  (2 
in.)  and  the  stand-offs  for  the  images  of  Figure  48  go  from  2.5  cm  (1  in.)  to  10.2  cm  (4  in.). 

Real  Part  Imaginary  Part 


Stand-off:  2.5  cm  5.1cm  7.6  cm  10.2  cm  2.5  cm  5.1cm  7.6  cm  10.2  cm 

Figure  48.  Effect  of  stand-off  distance  or  lift-off  on  array  response  for  vertical  scans  down  the  middle  of  cells  A3- 
A6.  The  vertical  axis  gives  the  scan  distance  in  inches.  The  horizontal  axis  gives  the  distance  along  the  sense 
element  array  and  spans  0.41  m  (16  in.).  The  object  type,  depth,  azimuth  angle,  dip  angle  and  approximate  position 
are:  M55  (0. 1  m,  0°,  45°)  at  ~28  in.,  M385  (0. 1  m,  0°,  45°)  at  -106  in.,  BLU-26  (0. 1  m,  0°,  0°)  at  -184  in.,  and  M86 
(0.4  m,0°,  45°)  at -262  in. 
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From  the  scan  images,  the  values  for  the  real  and  imaginary  parts  of  the  response  can  be 
extracted  for  the  center  of  the  object.  This  corresponds  to  the  object  location  directly  beneath  a 
sense  element.  These  values,  plotted  in  Figure  49,  show  how  the  response  approaches  the  “no 
object”  point,  which  is  the  origin  (0,0)  in  this  plot,  as  the  standoff  is  increased.  The  response  for 
the  BLU-26  is  approximately  linear,  similar  to  the  laboratory  measurements.  For  the  M385, 
which  tends  to  have  a  small  imaginary  response  at  this  measurement  frequency,  there  is  more 
sensitivity  to  noise  in  the  imaginary  part  of  the  signal.  Similarly,  for  the  M55,  which  tends  to 
have  a  small  real  response  at  this  measurement  frequency,  there  is  more  sensitivity  to  noise  in  the 
real  part  of  the  signal.  Furthermore,  as  will  be  shown  later,  the  sensor  proximity  to  the  ground 
tends  to  have  a  greater  effect  on  the  real  part  than  the  imaginary  part,  which  helps  to  explain  the 
negative  values  obtained  for  the  real  part  of  the  M55  response.  This  motivates  the  use  of  an 
additional  short  spatial  wavelength  sensor  in  the  follow-on  that  will  allow  the  local  elfects 
of  the  ground  (i.e.,  standoff  distance)  to  be  determined  so  that  it  can  be  compensated  for  in 
the  characteristic  response  of  the  object. 


Figure  49.  Plot  of  response  from  the  center  of  the  objects  using  the  scan  images  of  Figure  47. 

Measurement  scans  were  performed  over  a  number  of  the  other  grid  cells  as  well.  Magnitude 
images  resulting  from  horizontal  scans  across  columns  are  shown  in  Figure  50  and  Figure  51. 
Magnitude  images  resulting  from  vertical  scans  across  lanes  are  shown  in  Figure  55.  The 
emphasis  of  these  preliminary  scans  was  on  the  shallower  objects  to  demonstrate  the  imaging 
capability  of  the  array  and  to  determine  sensitivity  to  objects  using  a  prototype  drive  winding. 

From  the  horizontal  scans  of  Figure  50  and  Figure  51,  shallow  objects  tend  to  be  detectable  with 
this  array.  The  M385  and  BLU-26  are  evident  at  depths  of  0. 1  m  and  0.2  m  while  the  BDU-28  is 
apparent  at  the  0. 1  m  depth.  Although  UXO  deeper  than  0.2  m  tend  to  be  less  apparent  with  this 
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drive  construct,  some  of  the  larger  objects  are  apparent  at  deeper  depths.  For  example,  the 
M49A3  is  visible  at  a  depth  of  0.25  m,  as  indicated  by  Figure  53. 


D4  M385  (0.1  m,  n/a,  -90°)  E4  M385  (0.2  m,  0°,  0°)  F4  M385  (0.2  m,  90°,  0°)  G4  8#  shot  (0.2  m,  0°.  0°) 

Figure  50.  Horizontal  scans  over  the  middle  of  cells  in  Lane  4.  A  logarithmic  color  scale  for  the  magnitude  images 
is  shown.  The  depth,  azimuth  angle,  and  dip  angle  for  each  object  is  indicated.  The  horizontal  axis  gives  the  scan 
distance  in  inches.  The  vertical  axis  gives  the  distance  along  the  sense  element  array  and  spans  0.41  m  (16  in.). 


G5  8#  shot  (0.2  m,  0°,  0°)  H5  BDU-28  (0.1  m.  O',  45°)  15  BDU-28  (0.1  m,  0°, -45°)  J5  BDU-28  (0.1  m,  0°,90°) 


J5  BDU-28  (0.1  m,  0°,  90°)  K5  BDU-28  (0.1  m,  n/a, -90°)  L5  BDU-28  (0.2  m,  0°,  0°)  M5  BDU-28  (0.2  m,  90°,  0°) 

Figure  51.  Horizontal  scans  over  the  middle  of  cells  in  Lane  5.  A  logarithmic  color  scale  for  the  magnitude  images 
is  shown.  The  depth,  azimuth  angle,  and  dip  angle  for  each  object  is  indicated.  The  horizontal  axis  gives  the  scan 
distance  in  inches.  The  vertical  axis  gives  the  distance  along  the  sense  element  array  and  spans  0.41  m  (16  in.). 


Scans  were  also  performed  over  the  grid  cells  containing  clutter.  Figure  52  shows  the  results  of 
one  scan  across  the  center  region  of  the  cells,  containing  an  M49A3  shell  and  different  levels  of 
clutter.  The  intensity  of  the  measurement  signal  increases  with  the  amount  of  material  in  the  grid 
cell  and  is  greatest  for  the  cell  with  high  clutter.  The  color  scales  of  the  images  are  also 
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dominated  by  the  high  clutter  response  so  that  the  response  from  the  M49A3,  which  is  relatively 
deep  for  this  sensor  to  detect,  is  not  readily  visible. 


Figure  52.  Horizontal  scan  over  the  middle  of  cells  H8-K8.  All  include  M49A3  (0.25  m  depth,  0°  azimuth,  0°  dip), 
with  no  clutter,  low  clutter,  medium  clutter,  high  clutter  spaced  approximately  2  m  (78  in.)  apart  Top.  magnitude; 
Middle:  Real  part;  Bottom:  Imaginary  part  The  axes  give  the  scan  dimensions  in  inches.  The  vertical  axis  gives 
the  distance  along  the  sense  element  array  and  spans  0.41  m  (16  in.) 


Figure  53  shows  the  results  of  repeated  scans  over  the  areas  with  clutter  and  uses  a  logarithmic 
color  scale.  This  color  scale  lets  the  relatively  weak  signal  from  the  M49A3  itself  as  well  as  the 
strong  signals  from  the  clutter  to  be  visible.  The  color  scaling  also  accentuates  any  measurement 
noise  and  element-to-element  variation  in  the  array  to  be  more  noticeable.  The  scans  over  the 
cells  are  repeatable  but  also  show  some  drift  in  the  response  as  the  time  from  calibration  (e  g., 
after  repeated  scans  without  recalibration)  increases. 


M49A3  -  no  clutter  M49A3  with  low  clutter  M49A3  with  medium  clutter  M49A3  with  high  clutter 


Figure  53.  Magnitude  scans  only,  logarithmic  scale,  for  repeated  scans  over  the  grid  cells  with  clutter. 


Page  56  of  66 


The  calibration  area  containing  the  0.3  m  steel  plates  was  basically  underwater  and  inaccessible 
for  the  preliminary  measurements  with  this  system.  While  portions  of  the  area  around  the  0.6  m 
steel  plates  also  contained  water,  this  area  was  more  readily  accessible.  Figure  54  shows  the 
results  of  several  measurement  scans  over  the  0.6  m  steel  plates.  The  shallower  plate,  at  a  depth 
of  0.5  m,  is  readily  visible  and  provides  a  large  spatial  response.  Offset  scans  over  the  top  and 
bottom  section  of  the  grid  cells  capture  the  spatial  extent  of  the  plate  response  perpendicular  to 
the  scan  direction.  For  imaging  of  large  spatial  responses  such  as  this  a  wider  array  of  sensing 
elements  is  preferred. 
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Figure  54.  Horizontal  scans  over  cells  FI  1-H1 1 .  0.6  m  steel  plates  were  buried  to  a  depth  of  0.5  m  (near  a 
horizontal  position  of  70  inches)  and  to  a  depth  of  1.0  m  near  a  horizontal  position  of  236  inches.)  The  middle 
imagP!  was  a  scan  down  the  center  of  the  lane.  The  top  and  bottom  images  had  the  array  offset  by  0.3  m  (12  in.). 
The  axes  give  the  scan  dimensions  in  inches.  The  vertical  axis  gives  the  distance  along  the  sense  element  array  and 
spans  0.41  m  (16  in.) 


Vertical  scans  down  a  column  were  also  performed  over  some  of  the  objects  since  the  MWM- 
Array  is  directional  and  there  can  be  additional  information  gained  from  scanning  in  other 
orientations.  In  the  follow-on  efforts,  two  or  three  drive  orientations  will  be  used  in  an 
integrated  approach.  This  was  demonstrated  in  the  laboratory  measurements  showing  the 
orientation  of  the  several  objects.  The  results  of  some  of  these  scans  are  shown  in  Figure  55. 
Note  that  the  color  scale  was  chosen  for  these  figures  to  highlight  weak  responses  (e.g.,  M55  at 
0.1  m)  along  with  stronger  responses  (M385  at  0.1  m),  even  though  this  highlights  ground 
variations  and  noise  effects  as  well.  Figure  56  shows  the  real  and  imaginary  parts  of  the 
response  for  several  offset  scans  over  the  high  clutter  cell.  These  images  are  consistent  with  the 
horizontal  scans  of  Figure  52.  The  response  of  the  clutter  is  much  stronger  than  the  response  of 
the  single  objects,  but  also  spans  a  large  portion  of  the  grid  cell.  Additional  scans  were 
performed  with  the  array  offset  to  the  left  and  right  sides  as  well.  Even  though  these  offset  scans 
were  not  registered  very  well  with  respect  to  the  scan  down  the  middle  of  the  lanes,  they  provide 
an  indication  of  the  spatial  extent  of  the  clutter.  JENTEK  often  uses  spatial  signatures  to 
suppress  background  and  instrument  noise.  These  should  be  used  in  the  follow-on  effort 
forUXO. 
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Figure  55.  Vertical  scans  down  the  middle  of  several  grid  cells.  Object  type,  depth,  azimuth  angle,  and  dip  angle 
are  specified.  The  axes  give  the  dimensions  in  inches.  The  horizontal  axis  gives  the  distance  along  the  sense 
element  array  and  spans  0.41  m  (16  i.n). 
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Figure  56.  Results  of  three  vertical  scans  over  the  high  clutter  cell  (K8),  which  also  contains  an  M49A3  at  a  depth 
of  0.25  m,  0°  azimuth,  0°  dip.  The  scans  down  the  middle,  left,  and  right  sides  provide  spatial  information  about  the 
detected  object  over  a  wide  area. 


A  variety  of  baseline  or  diagnostic  measurements  were  also  performed.  As  an  example.  Figure 
57  shows  the  result  of  measurements  performed  with  a  stationary  array  over  a  grid  cell  that 
should  not  have  any  objects.  The  array  was  calibrated  prior  to  performing  the  measurement  of 
Figure  57a  and  the  plot  of  the  transient  data  versus  element  number  shows  that  any  drift  and 
element-to-element  variability  is  modest.  Figure  57b  shows  a  similar  plot,  without  a 
recalibration,  that  was  taken  after  a  scan  across  several  grid  cells.  While  the  drift,  indicated  by 
the  error  bars,  remains  modest,  there  is  noticeable  element-to-element  variability.  This  element- 
to-element  variability  is  attributed  to  the  mechanical  motion.  In  the  follow-on  effort,  methods 
to  reduce  drift  and  perform  rapid  recalibrations  should  be  included. 

Measurements  were  also  performed  over  grid  cells  that  did  not  contain  any  objects.  Figure  58 
shows  the  scan  data  for  one  set  of  grid  cells.  The  real  part,  and  hence  the  magnitude,  of  the 
response  shows  relatively  large  variations  near  the  center  area  of  the  scan.  In  contrast  the 
imaginary  part  does  not  show  significant  variation  along  the  scan;  the  striping  in  the  image 
reflects  the  element-to-element  variability  of  the  response.  The  response  near  the  center  area  is 
associated  with  variations  in  the  sensor  lift-off  from  the  ground  or  with  ground  property 
variations.  This  response  is  repeatable.  As  a  contrast,  a  scan  performed  over  another  group  of 
grid  cells  without  objects  is  shown  in  Figure  59  and  does  not  have  the  gradual  response  near  the 
center  of  the  scan.  This  scan  is  relatively  clear,  except  for  some  drift  in  the  real  part  of  the 
response  and  an  occasional  rapid  change  in  the  response  for  all  of  the  sense  elements.  This  type 
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of  response  can  be  removed  using  a  filter  that  removes  the  common  mode  response  at  a 
particular  measurement  location. 

For  reference.  Figure  60  shows  the  effect  of  lifting  the  array  over  a  fixed  grid  cell  location 
(without  an  object)  so  that  the  effects  of  proximity  to  the  ground  can  be  quantified.  This 
indicates  that  the  standoff  from  the  ground  can  have  a  significant  effect  on  the  real  part  of  the 
response,  but  a  minimal  effect  on  the  imaginary  part.  The  error  bars  in  the  plot  denote  the 
standard  deviation  of  the  measurements  at  each  height  and  across  all  of  the  sense  elements. 


Figure  57.  The  variation  in  response  between  dements,  where  the  error  bars  denote  the  standard  deviation  of  the 
measurements.  These  plots  reflect  no  motion,  stationary,  measurements,  (a)  The  array  was  recalibrated  prior  to 
these  measurements,  (b)  Data  was  taken  after  a  scan  over  approximately  four  grid  cells. 
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Figure  58.  Scan  over  grid  cells  without  objects,  in  the  middle  of  cells  between  Lanes  12  and  13  from  Column  F  to 
~J.  Top:  magnitude;  Middle:  Real  part;  Bottom:  Imaginary  part.  The  horizontal  axis  gives  the  scan  distance  in 
inches.  The  vertical  axis  gives  the  distance  along  the  sense  element  array  and  spans  16  inches  (0.41  m). 


Figure  59.  Scan  over  another  group  of  grid  cells  without  objects,  in  the  middle  of  cells  between  Lanes  2  and  3  from 
Column  between  E  and  F  to  ~J.  Done  without  rails.  Top:  magnitude;  Middle:  Real  part;  Bottom:  Imaginary  part 
The  horizontal  axis  gives  the  scan  distance  in  inches.  The  vertical  axis  gives  the  distance  along  the  sense  element 
array  and  spans  16  inches  (0.41  m). 
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Figure  60.  No  scan  motion  measurement.  The  array  was  raised  from  2.5  cm  (1  in.)  in  -15.2  cm  (6  in.)  increments 
every  500  points.  The  average  signal  is  from  all  16  elements  at  each  height 

Diagnostic  measurements  were  planned  for  the  test  pit,  where  scans  were  to  be  performed  over 
known  objects  at  specific  depths.  However,  the  conditions  of  the  test  pit,  shown  in  Figure  61, 
prevented  these  measurements  from  being  performed.  Water  was  drained  from  the  pit  on  the 
first  day  with  the  hope  that  it  would  dry  out  sufficiently  by  the  last  day  of  the  trial.  In  addition, 
significant  settling  and  dirt  transfer  outside  the  pit  boundaries,  where  digging  is  not  allowed, 
made  the  surface  of  the  test  pit  much  lower  than  the  surrounding  soil. 


Figure  61.  Photographs  of  the  test  pit  area. 
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Criteria  for  Enhancements  and  Recommendations 

This  base  program  was  aimed  at  generating  a  capability  baseline  for  the  development  of  an 
electromagnetic  induction  sensor  array  for  UXO  detection  and  clutter  suppression.  The  eventual 
goal  for  the  program  is  a  robust  discriminator  for  buried  objects.  Preliminary  measurements  in 
the  laboratory  and  the  field  were  performed  in  the  base  program  using  an  existing  prototype 
drive  winding,  with  minor  modifications,  and  an  array  of  sensing  elements  for  imaging  of  hidden 
objects.  Based  on  these  preliminary  measurements,  enhancements  to  the  system  have  been 
identified.  These  criteria  address  both  performance-based  and  instrumentation  issues,  as: 

•  Depth  of  sensitivity  to  UXO  —  This  reflects  the  signal-to-noise  ratio  of  the 
instrumentation  for  specific  objects  at  prescribed  depths.  The  noise  includes  background 
variations  in  soil  properties  as  well  as  drift  and  instrumentation  noise.  Another  aspect  of 
this  criterion  is  the  accuracy  and  robustness  of  the  estimates  of  the  object  properties. 

•  Multiple  object  resolution  —  This  reflects  the  capability  to  isolate  or  characterize  the 
response  of  individual  objects  that  are  near  one  another.  This  could  include  benign 
clutter  over  a  UXO  or  multiple  proximate  UXO.  Similar  to  the  depth  of  sensitivity,  this 
could  be  defined  in  terms  of  specific  objects  at  prescribed  depths  and  with  prescribed 
separations. 

•  Sensitivity  to  soil  properties  -  Variations  in  soil  properties,  such  as  the  magnetic 
permeability  or  electrical  conductivity,  can  affect  the  response  of  electromagnetic 
induction  instruments  and  can  mask  the  response  of  UXO  unless  accounted  for. 

•  Probability  of  detection  -  This  provides  a  measure  of  system  performance  for  detecting 
all  relevant  objects. 

•  Confidence  of  identification  -  This  addresses  the  capability  to  identify  an  object  as 
harmful  or  benign  after  detection. 

•  Data  acquisition  rate  and  scan  speed  —  Efficient  site  reclamation  requires  sufficiently 
high  data  acquisition  and  scan  rates  for  the  signature  responses  of  the  buried  objects  to  be 
captured  and  characterized. 

•  Instrument  bandwidth  —  This  should  be  large  enough  so  that  measurements  capture 
enough  of  the  frequency  dependent  response  to  assist  in  the  characterization  of  the  UXO. 
Higher  bandwidth  may  reduce  the  achievable  scan  speed. 

•  Instrument  portability  -  The  size,  weight,  and  cost  of  the  instrument  are  considerations 
for  deployment  and  field  use. 

•  Sensor  array  scalability  -  The  array  needs  to  be  large  enough  to  provide  sensitivity  to 
deep  objects  within  the  typical  depth  of  UXO  (e.g.,  1  m). 

•  Data  processing  —  Data  obtained  for  the  objects  should  be  processed  rapidly  and 
accurately  to  provide  an  assessment  of  the  object  characteristics.  This  data  includes,  for 
example,  the  scan  images  as  well  as  multiple  frequency,  sense  element  location  and  sense 
element  orientation  information. 

•  Observability  of  object  features  -  the  ability  to  independently  measure  multiple  object 
features  such  as  size,  depth,  orientation,  material  type,  and  shape  (e.g.,  aspect  ratio). 
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Recommendations  for  the  follow-on  effort  include: 


•  Deeper  penetration  depth  drive  -  The  drive  winding  should  be  modified  to  permit  a 
deeper  penetration  depth  for  the  magnetic  field  into  the  soil.  The  prototype  drive  used  in 
this  base  program  has  a  depth  limitation  of  roughly  0.45  m  (18  in.),  including  the  standoff 
distance  to  the  ground  surface,  so  the  focus  has  been  on  objects  located  within 
approximately  0.4  m  (16  in.)  of  the  surface.  The  depth  for  UXO  can  be  up  to  1  m. 

•  Wider  drive  and  larger  array  of  sense  elements  -  Since  edge  effects  were  visible  at  the 
outermost  sense  elements  of  the  array,  to  improve  the  image  quality  the  drive  winding 
should  be  modified  by  extending  the  drive  winding  segments  further  past  the  sense 
element  array.  A  longer  array  of  sense  elements  facilitates  imaging  of  deeper  and  larger 
objects  since  a  single  scan  of  the  array  may  provide  a  complete  spatial  image  of  the 
response.  A  higher  resolution  array  (e.g.,  32  instead  of  16  elements)  would  also  improve 
the  capability  to  determine  object  shape. 

•  Multiple  spatial  wavelength  measurements  -  Modifying  the  sensor  design  to  perform 
measurements  at  multiple  spatial  wavelengths  would  permit  compensation  for  local 
ground  property  and  stand-off  variations.  This  should  also  permit  the  separation  of  near 
surface  object  and/or  soil  property  effects  on  the  response  from  deep  object  property 
variations.  JENTEK  often  uses  spatial  signatures  to  suppress  background  and  instrument 

noise. 

•  Addition  of  complementary  sense  elements  -  Sensing  elements  placed  at  multiple 
locations  within  the  drive  footprint  and  at  different  orientations  capture  complementary 
information  about  the  buried  object.  This  could  be  accomplished  with  a  two-dimensional 
array  of  sense  elements. 

•  Reduced  sense  element  variability  -  This  is  necessary  for  the  construction  of  high 
resolution  images  and  can  be  accomplished  with  hardware  and  software  modifications. 
For  example,  drift  can  be  reduced  with  occasional  rapid  recalibration.  In  software, 
normalizing  across  the  elements  at  a  user-selected  location,  preferably  not  over  an  object, 
would  reduce  image  striping. 

•  Improved  ease  of  use  -  The  system  should  be  more  portable  for  field  measurements  and 
may  include  mounting  equipment  onto  the  cart. 

•  Reduced  mechanical  sensitivity  of  array  to  cart  motion  -  The  cart  design  should  be 
modified  to  reduce  tilting  of  the  array  over  uneven  terrain  while  maintaining  the  ability  to 
raise  and  lower  the  array  over  detected  objects.  Filters  can  remove  the  common  mode 
(noise)  response  at  a  particular  measurement  location. 

•  Additional  location  information  -  Incorporation  of  a  GPS  locator  or  some  other 
location  method  would  facilitate  scanning  of  large  areas.  The  use  of  a  single  position 
encoder  can  provide  accurate  local  position  measurements  but  only  in  one  direction. 

•  Multiple  frequency  measurements  -  The  spectral  response  reflects  material  and 
geometric  properties  of  the  UXO.  Measurements  have  been  performed  at  a  single 
frequency,  but  there  is  additional  information  that  could  be  obtained  from  multiple 
frequency  measurements. 

•  Detection  and  suppression  algorithm  enhancements  -  The  algorithms  for 
characterizing  the  objects  and  clutter  suppression  need  to  integrate  the  responses  from 
multiple  sensing  elements,  preferably  at  multiple  element  locations  and  orientations 
within  the  drive  footprint,  with  multiple  measurement  frequencies.  They  also  need  to 
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combine  information  from  high  resolution  images  with  property  estimates  from 
measurement  grids  or  some  other  model  based  approach.  Measurement  responses  from 
two  or  three  drive  orientations  should  also  be  incorporated. 

•  Nonspherical  object  model  development  -  Basic  object  property  estimates  have  been 
obtained  with  a  sphere-based  model,  but  better  estimates  are  expected  when  the 
nonspherical  shape  of  typical  UXO  is  accounted  for.  The  use  of  ellipsoid  models  with 
multiple  sensing  elements  to  discriminate  object  features  (shape,  size,  orientation)  should 
be  investigated. 

•  Multiple  drive  orientations  -  Multiple  drive  winding  orientations  are  necessary  to 
improve  observability  for  different  object  types. 
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Concluding  Remarks 


This  base  program  was  aimed  at  generating  a  base-line  response  for  an  electromagnetic  induction 
imaging  array  for  UXO  detection  and  discrimination.  The  focus  was  on  adapting  our  impedance 
instrumentation  to  function  with  a  prototype  array  to  demonstrate  capabilities  and  limitations  of 
this  early  stage  system  and  to  demonstrate  the  use  of  models  for  extracting  information  about 
object  features.  All  technical  milestones  for  this  base  program  have  been  met.  These  milestones 
led  to  the  generation  of  enhancement  criteria  and  recommendations  that  are  expected  to  be 
pursue  in  the  follow-on  program. 

The  base  program  accomplishments  include  (1)  adapting  an  existing  drive  winding  and  an  array 
of  sense  elements  to  work  with  parallel  architecture  impedance  instrumentation;  (2) 
demonstrating  that  improving  (e.g.,  doubling)  the  array  resolution  can  be  accomplished  without 
affecting  the  signal-to-noise  ratio;  (3)  demonstrating  a  rapid-scanning  capability;  (4) 
demonstrating  that  scan  images  can  reflect  object  orientation  and  shape  and  also  the  presence  of 
multiple  proximate  objects;  (5)  demonstrating  a  preliminary  size  and  depth  estimation  capability 
using  a  model-based  approach;  and  (6)  a  similarity  of  object  signature  responses  in  the  laboratory 

and  the  field. 

System  enhancements  recommended  for  the  follow-on  effort  include:  (1)  modifying  the  drive 
winding  for  deeper  object  sensitivity  and  to  reduce  edge  effects  on  the  image;  (2)  fabrication  of  a 
two-dimensional  array  of  sense  elements;  (3)  multiple  spatial  wavelength  measurements,  (4) 
algorithm  enhancements  to  suppress  background  effects,  sense  element  variability,  and  clutter; 
(5)  methods  to  reduce  common  mode  noise;  and  (6)  use  of  two  or  more  drive  winding  and 
sensing  element  orientations  to  further  improve  observability. 
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